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Abstract 

In this paper we make a further step towards a dispersive description of the hadronic light-by-light (HLbL) tensor, 
which should ultimately lead to a data-driven evaluation of its contribution to [g — 2)^. We first provide a Lorentz 
decomposition of the HLbL tensor performed according to the general recipe by Bardeen, Tung, and Tarrach, generalizing 
and extending our previous approach, which was constructed in terms of a basis of helicity amplitudes. Such a tensor 
decomposition has several advantages: the role of gauge invariance and crossing symmetry becomes fully transparent; the 
scalar coefficient functions are free of kinematic singularities and zeros, and thus fuffill a Mandelstam double-dispersive 
representation; and the explicit relation for the HLbL contribution to [g — 2)^ in terms of the coefficient functions 
simplifies substantially. We demonstrate explicitly that the dispersive approach defines both the pion-pole and the 
pion-loop contribution unambiguously and in a model-independent way. The pion loop, dispersively defined as pion-box 
topology, is proven to coincide exactly with the one-loop scalar QED amplitude, multiplied by the appropriate pion 
vector form factors. 
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1 Introduction 

The anomalous magnetic moment of the muon, = (g — 2)^/2, is one of the very rare quantities in particle 
physics where a significant discrepancy between its experimental determination and the Standard-Model eval¬ 
uation still persists. A robust, reliable estimate of the theoretical uncertainties is thus mandatory, especially in 
view of the projected accuracy of the forthcoming {g — 2)^ experiments at FNAL [1] and J-PARC [2], to settle 
whether this discrepancy can be attributed to New Physics. 

The uncertainty on the theory side is dominated by hadronic contributions, see e.g. [3-5], the largest error 
arising from hadronic vacuum polarization (HVP). Since HVP is straightforwardly related to the total e“'"e“ 
hadronic cross section via a dispersion relation, the evaluation of this contribution is expected to become more 
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accurate within the next few years [ 6 ] thanks to upcoming improvements in the experimental input/ although 
reducing further the present sub-percent accuracy will be challenging. This implies that hadronic light-by-light 
(HLbL) scattering is likely to soon dominate the theory uncertainty in Using previous approaches [11- 
22], systematic errors are difficult, if not impossible, to quantify. A novel strategy is required to go beyond 
the state of the art, to avoid or at least reduce model dependence to a minimum, make a solid estimate of 
theoretical errors, and possibly reduce them. Lattice QCD is a natural candidate to achieve this goal, but 
it is not clear yet if or when this method will become competitive [23-25]. Alternatively, as we recently 
argued [26], it is possible to establish a rigorous framework based on dispersion theory that directly links 
to experimentally accessible on-shell form factors and scattering amplitudes [27], contrary to what has been 
previously conjectured [12, 24, 28, 29].^ 

Our dispersive formalism is based on the fundamental principles of unitarity, analyticity, crossing symmetry, 
and gauge invariance. The derivation of such a systematic framework becomes challenging due to the fact 
that HLbL scattering is described by a hadronic four-point function whose properties are significantly more 
complicated than those of the two-point function entering HVP. Besides the opportunity to achieve a data- 
driven evaluation of , this approach allows us to unambiguously define and evaluate the various low-energy 
contributions to HLbL scattering, most notably pion pole and pion loop. 

The key result of [26] was a master formula giving the contribution to from the pion-pole term and 

intermediate tttt states expressed in terms of the pion transition form factor and helicity partial waves for 
7 * 7 * —> TTTT. In the present paper we reformulate this approach in a more general context, by employing a 
different generating set for the Lorentz structures of the HLbL tensor. This new set, which we constructed 
following the prescription by Bardeen, Tung [31], and Tarrach [32] (BTT), has the property of being explicitly 
free of kinematic singularities and zeros [33]. Therefore, the scalar coefficient functions associated with each 
Lorentz structure of this generating set are also free of kinematic singularities and zeros, and thus fulfill a 
Mandelstam double-dispersive representation. Further motivation for deriving such a BTT set is provided by 
the fact that the role of gauge invariance and crossing symmetry becomes fully transparent, with constraints 
from soft-photon zeros incorporated automatically. Moreover, the absence of kinematic singularities makes it 
much easier to perform the angular averages necessary to calculate the contribution to {g — 2 )^ — indeed both 
the derivation and the final form of the master formula are simplified considerably in the present setting, and 
allow for a simpler inclusion of partial-wave contributions beyond 5'-waves. In the derivation of the master 
formula it is very useful to perform a Wick rotation of the integration variables. This has been already used in 
several papers in the literature without ever addressing an important subtlety which is relevant in the present 
case: namely if the Wick rotation goes through without changes even if the integrand is a loop function showing 
anomalous thresholds. We discuss this point in detail and prove that this is indeed the case. Finally, the BTT 
formulation facilitates the proof that both pion-pole and pion-loop contributions are uniquely defined in the 
dispersive approach, and confirms the relation to the pion transition form factor and the pion electromagnetic 
form factor given in [26]. We explicitly proof that the scalar QED (sQED) pion loop dressed with pion form 
factors (denoted by FsQED in [26]) is identical to the contribution from tttt intermediate states with a pion-pole 
left-hand cut (LHC). 

In summary, this paper contains three main results: (i) the explicit construction of the BTT decomposition 
of the HLbL tensor; (ii) the derivation of the master formula expressing the HLbL contribution in terms of 
the BTT scalar functions; (iii) the proof that in the dispersive approach the FsQED contribution is uniquely 
defined. Everything which goes beyond the FsQED contribution is amenable to partial-wave expansion and will 
be discussed in full detail in a forthcoming publication. 

The present paper is structured as follows. In Sect. 2, we first review some aspects of the process 7 * 7 * —t 
TTTT, mainly as an illustration of the techniques that we apply afterwards to HLbL. In Sect. 3, we derive the 
decomposition of the HLbL tensor into a set of scalar functions that are free of kinematics. In Sect. 4, this BTT 
decomposition is then used to derive a master formula for that is parametrized by the corresponding 

scalar functions. We verify the validity of the Wick rotation of the two-loop integral in the calculation of 
even in the case of anomalous thresholds in the HLbL tensor. In Sect. 5, we derive a Mandelstam representation 
of the BTT scalar functions, and, by comparison to the loop formulation of the sQED pion loop, demonstrate 
that it indeed coincides with the contribution of tttt intermediate states with a pion-pole LHC. We conclude with 
an outlook in Sect. 6 , while some of the lengthier expressions and derivations are collected in the appendices. 

^For a recently suggested alternative approach to determine the leading hadronic correction to a^, based on space-like data 
extracted from Bhabha scattering, see [7]. 

^At this order in the fine-structure constant also two-loop diagrams with HVP insertion appear [8]. Higher-order hadronic 
contributions have been investigated in [9, 10]. 

^For a different approach aiming at a data-driven evaluation of based on a dispersive description of the Pauli form factor 

instead of the HLbL tensor, see [30]. 
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2 The sub-process 7 * 7 * —)■ tttt 

As a prelude, we discuss in this section the process 7 * 7 * —>■ tttt, which will become important as a sub-process 
when we write down a dispersion relation for the HLbL tensor. Because the Lorentz structure of 7 * 7 * —> tttt 
is much simpler than the one of light-by-light scattering, it also allows us to illustrate a technique for the 
construction of amplitudes that are free of kinematic singularities and zeros, which we will apply afterwards to 
the more complicated case of HLbL. In order to render this analogy manifest, we slightly change conventions 
compared to [26]: we identify here 7 * 7 * —>■ tttt as the s-channel process instead of 7 * 7 r —)■ 7 * 7 r. This assignment is 
slightly unnatural when it comes to constructing dispersion relations for a process with these crossing properties, 
see [34-37], but makes the inclusion into the unitarity relation in HLbL more straightforward. 

2.1 Kinematics and matrix element 



k2 ki 

Figure 1: 7 * 7 * —>■ tttt as a sub-process of e''"e“ —>■ e''"e“7r7r. 

Consider the process 

e+(fci)e“(A:2) e+(A:3)e“(A:4)7*(9i)7*(92) e+(fc3)e“(fc4)7r“(pi)7r*'(p2), ( 2 . 1 ) 

shown in Fig. 1. At O(e^), the amplitude for this process is given by 
iT = v(ki)(-ieja)v(k3)u(k4)(-ie^0)u(k2) 

X ^ ( 2 . 2 ) 

9i v 9i / 92 \ 92 y 

where ^ is an arbitrary gauge parameter for the photon propagators and the tensor is defined as the pure 

QCD matrix element 

W^biPi^P2,qi)='i J d'^xe-^‘^^'^{TT‘^{pi)TT\p2)\T{j^^{x)j^^{0)}\0). ( 2 . 3 ) 

The contraction thereof with appropriate polarization vectors may be understood as an amplitude for the 
off-shell process 

7*(9i, Ai)7*((72, A 2 ) 7r“(pi)7r*'(p2), (2.4) 

where Ai 2 denote the helicities of the off-shell photons. We define the connected part of such a matrix element 

by 

(7r“(pi)7r'’(p2)|7*(9i, Ai)7*(92, A 2 )) 

:= - eh^^Hqi)4Hq2) J Ae-*(^-"+^-")(7r7pi)7r'’(p2)|r{je^,n(^)jem(y)}l^^ 

= - e727r)^(5('‘)(pi -bp 2 - 9i - 92)e^K9i)e^" ( 92 ) 

X J d'^xe-^‘>^'^{7r‘^{pi)7r\p2)\T{j^^{x)j^^{0)}\0) 

=ie'^{2Tr)‘^S^‘^'>{pi + P 2 - qi - {pi,p 2 ,qi). ( 2 . 5 ) 
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The helicity amplitudes are given by the contraction with polarization vectors: 


We introduce the following kinematic variables:"^ 

s := (<71 +q 2 f = (pi +P 2 )^ 
t := {qi -pif = {q 2 -p 2 f, 
u := (<71 -P2f = {q2-Pif, 


( 2 . 6 ) 


(2.7) 


which satisfy s + t + u = ql+q 2 + 2Mj. 

2.2 Tensor decomposition 

The tensor can be decomposed based on Lorentz covariance as (we drop isospin indices) 

+ gf (2.8) 

where we abbreviate qi = {qi,q 2 ,P 2 — Pi} and where double indices are summed. The ten coefficient functions 
{Wi, cannot contain any kinematic but only dynamic singularities. However, they have to fulfill kinematic 
constraints that are required e.g. by gauge invariance, hence they contain kinematic zeros. Conservation of the 
electromagnetic current implies the Ward identities 

= <72 = 0, (2.9) 

hence a priori six relations between the scalar functions {Wi, W}^}. Only five of them are linearly independent 
and thus reduce the set of scalar functions to five independent ones. 

We now construct a set of scalar functions which are free of both kinematic singularities and zeros. We 
follow the recipe given by Bardeen, Tung [31], and Tarrach [32]. As was shown in [32], the basis (consisting 
of five functions) constructed according to the recipe of [31] is not free of kinematic singularities. However, a 
redundant set of six structures can be constructed, which fulfills the requirement. 

In the case of 7 * 7 * —>■ tttt, there exists in addition to gauge invariance and crossing symmetry of the photons 
also the crossing (Bose) symmetry of the pions, which is responsible for additional kinematic zeros. In fact, this 
can be used to circumvent the introduction of a sixth function, reducing the set again to five scalar functions 
free of both kinematic singularities and zeros. This was first derived in [38] in the context of doubly-virtual 
nucleon Compton scattering.^ 

We define the projector 


JMXx (2.10) 

qi ■ 92 

which satisfies 

= W^xl\ = ( 2 . 11 ) 

i.e. the tensor is invariant under contraction with the projector, but contracting the projector with any 
Lorentz structure produces a gauge-invariant structure. 

We apply this projector for both photons: 

5 

( 2 . 12 ) 

“^In this Sect. 2 only, s, t, and u refer to the Mandelstam variables of 7 * 7 * tttt and not to the ones of the HLbL tensor. 
^This does not mean that the Lorentz structures are independent in all kinematic limits, hence Tarrach’s conclusion that no 
minimal basis exists is not invalidated. We thank J. Gasser for bringing this to our attention. 
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where 



M y 

92 9i 

9i - 92 ’ 


Bi = 

= 9i92 

9? 92 92 

9i9i9i 1 9?9i92 9i 

B2 = 

9i - 92 

9i - 92 (91 - 92)2 ’ 

= 9i93 

9? 92 93 

9i - 92 

92 - 939^91 1 9i92-939291 

9i - 92 (91 - 92)2 

B3 = 

II 

00 

to 

9^93 91 

91-9392 92 1 9i9i-939291 

Ba = 

9i - 92 

9i - 92 (91 - 92)2 

Li y 

= 9393 

9i - 9392 93 92 - 9393 91 . 9i - 9392 - 9392 91 

B5 = 

9i - 92 

9i - 92 (91 - 92)2 


rl3 


/32 
^2 ’ 


/33 


(2.13) 


As the functions Bi are a subset of the original scalar functions, they are still free of kinematic singularities, 
but contain zeros, because the Lorentz structures contain singularities. We have to remove now the single and 
double poles in qi ■ q 2 from the Lorentz structures . This is achieved as follows [31]: 

• remove as many double poles as possible by adding to the structures linear combinations of other structures 
with non-singular coefficients. 


• if no more double poles can be removed in this way, multiply the structures that still contain double poles 
by qi ■ (? 2 , 


• proceed in the same way with single poles. 

It turns out that no double poles in qi ■ q 2 can be removed by adding to the structures multiples of the other 
structures, hence , ■ ■ ■, have to be multiplied by qi ■ q 2 . The resulting simple poles can be removed by 
adding multiples of In the end, we have to multiply by qi ■ q 2 in order to remove the last pole. We 
then arrive at the following representation: 


where 


and 


W„ 


= Vf* 




Tr = qi ■ q 2 g^‘' - 92 91 > 

+ 9i • 929 i 92 - 9?92 92 “ 9^91 9i, 

^3 =9i 92 - 939 + 9i-929193 - 9i92 93 - 92-939191, 

Tr = 9291 - 939^'" + 9i - 9293 92 “ 9293 9i - 9i - 9392 92, 

= 9i - 9392 - 939^“" + 9i - 9293 93 “ 9i - 9392 93 “ 92 - 9393 9^, 


Bi 

B2 


Ba 


B5 


1 


9i 

- 92 


1 

9i 

- 92 


1 

9i 

- 92 


1 

9i 

- 92 


1 

9i 

- 92 


Wi - 

W^ 2 ^ 

W2^^ 

wf, 

wp. 


9?9i 

(91 - 92)^ 




12 


9192 - 93 

(91 - 92)^ 


wp 


qhl ■ 93 t^^ 32 _ 9l - 9392 - 93 tt/33 

( 91 - 92)2 " ( 91 - 92)2 


(2.14) 


(2.15) 


(2.16) 


The Lorentz structures (2.15) agree with the basis used in [26], upon q 2 —>■ — 92 , Pi —t —pi due to crossing and 
the identification 




_ /t-i[26J rp _ /t-i[26J 
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As shown by Tarrach, this basis is not free of kinematic singularities and zeros [32]. The structures 
form a basis for qi ■ q 2 ^ Q, but are degenerate for gi • 92 = 0. This degeneracy implies that there is a linear 
combination of the structures that is proportional to qi ■ q 2 ' 

qi ■ 9392 • 93Tr - 92^1 • - 9?92 • q.f^ + q^f^^ = q, ■ q2fr, (2.18) 


where 


Tr = {qlqz - 9i • 939 i ) ( 9I93 - 92 • 9392) • (2-19) 

Hence, in order to have a generating set even at gi • 52 = 0 this sixth structure has to be added by hand. 
Projected on the basis, it gives coefficients with poles in qi • ( 72 - Although the basis is not free of kinematic 
singularities and zeros, we have found the exact form of the singularities: 


5 6 

i=l i=l 

where 


( 2 . 20 ) 


Bi — Bi, 

5 D , 91 • 9392 • 93 

i>2 — i>2 H-i>6) 

9i • 92 

6 O 9 I 91 • 93 „ 

i>3 — i>3--D 6 , 

9i • 92 

6 o 9?92 • 93 „ 

r>4 — n4 -±> 6 , 


B,=B, + 


9i • 92 
2 2 
9192 ^ 

9i • 92' 


( 2 . 21 ) 


The functions are free of kinematic singularities. The tensor structures can be written in terms of the 
Mandelstam variables t and u as 


Tr = 9i • 929 '“' - 92 9i, 

Tr = qlqlq^’' + 9i • 929192 - 9?9292 - qlqiql, 

Tr = 91 • 929^93 - 9?9293 “ “ w)9?9'"" + ^(^ “ w)9i9r, 

'Tl'" = 9i • 9293 92 - 9I9391 + ^(^ - u)ql9^'' '“)92 92: 

Tr = 91 • 929393 - + ^(^ - w) (9391 - 92 9^ . 

= qlqlqaql + \it-u) (9^9392 -929193) - - ufqU2- (2.22) 

The derivation of the Lorentz decomposition for 7 * 7 * —> tttt up to this point should illustrate the techniques 
that will be applied to the much more complicated case of HLbL in Sect. 3. In the case of 7 * 7 * —>■ tttt, there is in 
fact a solution to get rid of the additional redundant structure without introducing kinematic singularities [38]. 
The necessary steps are discussed in the next subsection. Note that it is an open question whether an analogous 
procedure to get rid of all redundancies in HLbL exists, see the discussion below (3.29). 


2.3 Kinematic zeros due to crossing antisymmetry 

We have now constructed a redundant set of six Lorentz structures. All kinematic constraints from gauge 
invariance are implemented and the scalar coefficient functions Bi are free of kinematic singularities. However, 
due to the crossing symmetry of the pions, additional kinematic zeros are present in the functions Bi. As noted 
in [38], this can be exploited to eliminate the redundancy and to work again with basis coefficient functions free 
of kinematic singularities and zeros. 

Since the two-photon state is even under charge conjugation, so must be the two-pion state. Therefore, the 
isospin / = 1 amplitude vanishes. Bose symmetry implies that the amplitude and hence the tensor is 
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invariant under pi -f-)- p 2 or equivalently -f-)- —q^. Under this transformation, the tensor structures are 
even for i = 1,2,5, 6 and odd for i = 3,4. Hence, Hi 2 , 5,6 Eire even, while B 3 4 must be odd. They contain a 
kinematic zero of the form 


^3,4 = {t — u)B3^4, 


(2.23) 


where 53,4 are free of kinematic singularities. 

Crossing symmetry of the photons requires the invariance of under qi ^ q 2 , p ^ v. While T ^2 5 6 

invariant under this transformation, we observe the crossing relation -f-)- . Hence, for fixed Mandelstam 

variables. Hi, 2 . 5,6 and H 3 — H 4 are even under qi q^j while H 3 + H 4 is odd and contains a kinematic zero 
2 2 
qt - qi- 

We make a basis change 


where 




5 6 


pifliy 

pifliy 

pifliy 


_ p]^y 

_ pfiy 


= {t-u){f!i'' - fr ), 

_ pfiy 

_ rpiiy 


(2.24) 


(2.25) 


i.e. we trade off the combination , which is even under photon crossing, against the Tarrach structure 

Tq'^ and absorb the kinematic zero in t — m into . The projected basis functions are then: 


Ai — Hi, 

(t — u)^ H3 + H4 

2 ql - ql ’ 

B3 — B4 qi + 92 B3 + B4 


A 2 — H 2 + 
^3 = 


2 9? - 9| ’ 


H 4 = H 5 - 2 qlql ^l^^^ , 

qt - qi 

H5 = He + (s - 9? - ql) H^^2 ■ 

9i 92 


(2.26) 


The apparent kinematic singularities in q\ — are canceled by the corresponding zero in H 3 + H 4 . Hence, the 
basis functions Ai are free of both kinematic singularities and zeros. 


2.4 Helicity amplitudes and soft-photon zeros 

In the following, we construct the helicity amplitudes with the momenta and polarization vectors in the s-channel 
center-of-mass frame. We define the particle momenta as 


where 


9i = (Hqi,0,0, 19 I), 92 = (H,,,0,0,-|9|), 

Pi = (Ep, Ip I sin 0 cost/), |p| sinflsint/), |p| cos0), 

P2 = (Hp, —IpI sin0cos0, —IpI sin^sint/), —IpI COS0), 


Bqi = yqt + q^ = 


I 2 2 

S + 9i - 92 


) Bq2 = yqi + P = 


s-qt + qi 


Hp= VM2+p2 = ^, |p| = y^-M2 = ^a.(s). 


I9l = 


A^/ ^(5,9^9i) 

2 v^ 


(2.27) 


( 2 . 28 ) 




We use the notation 



The scattering angle is then given by 

2 := COS 0 

We define the polarization vectors by 


A(a, b, c) := + b^ + — 2{ab + bc + ca). 


t — u t — u 

4|g||p| CT^(s)Al/ 2 (s,g 2 ^g 2 )- 


e±(gi) = T;^( 0 , l,±i, 0 ), 
eo(<?i) = ^i\q\,0 ,0,EqJ, 
e±(<? 2 ) = T;^( 0 ,l,Tb 0 ), 

eofe) = ■^(-|g|,0,0,£;,J 

b .2 


(2.29) 


(2.30) 


(2.31) 


The longitudinal polarization vectors are normalized to one for = \/^. However, as the off-shell photons are 
not physical states, the choice of has no influence on any physical observable. The fact that the dependence 
on the normalization has to drop out in the end provides a useful check on the calculation, so that we will keep 
general in the following. 

The helicity amplitudes are then given by 


H++ = i?— = -i^(s - 91 - 92)211 - 9^92^2 + ;^(s - 4:M^)Xi2{s)z'^{qf + ql)A3 
2 2s 

+ i(s - 4M^) ((s - 9 ? - 9l) + - {ql + ,1)) A, 

+ - 4M2)(1 - z^)A3, 

H+- = H_+ = -^(s - 4M2)(1 - - 9i - 92)214 -t 2qlqlA3"j , 

H+o = -H_o = “ 4M^)zV 1 - 2 :^ (^Ai2(s)H3 - {s + qf - qDA^ - 9i(s - 9i + 91)^5^ , 

i?o+ = —Hq- = ~ 4 M,^)z \/1 — ^Ai2(s)H3 — (s — -|- q^)Ai — 92(5 + 9 i ~ 92)^5^ j 

Hoo = ( — Ai — -{s — ql — 92)212 -(s — 4 M^)Ai2(s)2:^H3 

C 16 V 2 s 

+ {s- 4M2)z^H4 + - 4M2) (s 2 _ (q ,2 _ ^ 2 ^ 2 ^ z^Hg^ , 


where Ai 2 (s) := A(s, 94 , ql). Since the functions Ai are free of kinematic singularities and zeros, we can read off 
from these equations the soft-photon zeros [39, 40]. In the limit 91 —> 0, the Mandelstam variables become 

s = ql, t = u = M^. (2.33) 


We conclude that the helicity amplitudes must vanish at this point apart from terms containing a dynamic 
singularity, which will be discussed in Sect. 2.6. The second soft-photon limit, 92 —t 0, leads to 

s = ql, t = u = Ml (2.34) 


and the same arguments apply for the helicity amplitudes. 

Crossing symmetry of the photons implies that under the transformation ql ql (and fixed Mandelstam 
variables), H++, _, and Hqq remain invariant, the other two helicity amplitudes transform as H+o ^ Ho+. 
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2.5 Partial-wave expansion 

To make the connection to [26], we briefly discuss the consequences of the decomposition (2.32) for the kernel 
functions that appear in partial-wave dispersion relations for 7 * 7 * —> tttt. Partial waves are most conveniently 
defined for helicity amplitudes. Using the formalism of [41], we write the partial-wave expansions as 


HaiAs (s, t, u) = J^( 2 J + l)dio(z)hj,A,x^ (s), (2.35) 

J 

where is the Wigner d-function, to = |Ai — A 2 I, and the helicity partial waves hj^x^x^ depend implicitly on 
the photon virtualities q\ and q^- Since the isospin of the two-pion system is / = 0, 2, only even partial waves 
are allowed. For to = 0, i.e. id++ and idoO) the partial-wave expansion starts at J = 0, otherwise at J = 2. 
Note also that doo(-z) = Pj{z) are the Legendre polynomials. 

The partial waves can be obtained by projection: 

1 

hj,AiA2(s) = J ^dzd:^Q{z)Hx^x2is,z). (2.36) 


Since the functions Ai are free of kinematic singularities and zeros, they are appropriate for a dispersive 
description. To this end, the Roy-Steiner treatment of [36] can be generalized to the doubly-virtual case. One 
starts by writing down hyperbolic dispersion relations 


A,{s,t,u) = ^f°™(s,t,'u) 


1 

JAMl s' — S 


1 /■“ 

-I — / dt'lmAi{t\u) 

^ Jto 


t' -t 


1 

t' — U 


t' — a 


(2.37) 


with hyperbola parameter a [34]. If we invert (2.32) to express the Ai in terms of the helicity amplitudes and 
insert the partial-wave expansion both on the left- and right-hand side of the dispersion relation, we obtain a 
set of Roy-Steiner equations 


hj,iis)=J2J2 


I AMI 


ds'Kyj,{s, s')lmhj'j{s') -|- ..., 


(2.38) 


where i,j G {1,2, 3,4, 5}, 


hj,i ■— hj ^2 ■— 

hj,3 ■= %hj+o -f %hjo+, hjA ■= 

qi ’ qf 


^2 1 7 7 _ _ ^ 1^2 

~~2 ^ J ,-\-0 ~~2 ^ ^ J ,5 •— ~2 2 

<?2 it qm 


2 2 hj,00j 


(2.39) 


and K'jj, are integral kernels. The ellipses in (2.38) stand for the contribution of the partial waves of the crossed 
channels. The normalization of the helicity amplitudes is reabsorbed in (2.39), so that the partial waves hj^i 
remain finite in the limit qf ^ 0. 

If only S'-waves are taken into account, the relation between the scalar functions and the helicity partial 
waves in the s-channel is given by 

2 

Ai = , . {2qlqlho,5{s) -{s-qf- ql)ho,i{s)) , 

Ai 2 (s) 

^2 = 7 , , (2ho,i(s) - (s - - g2)V5(s)) , 

Ai 2 (s) 

A3 = A4 = A5= 0. (2.40) 


In this case, the scalar functions depend only on s. If we take into account U-waves as well, the scalar functions 
A3, Aa, and A3 no longer vanish, but they depend only on s, while the scalar functions Ai and A2 are now 
second order polynomials in t and u. The explicit expressions are given in App. A. 
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For the diagonal kernel functions in (2.38), we find: 


K],\{s, s ) 
Kll{s,s') 

^22 (Sj *0 

K^2{s,s') 


Kll{s,s') 

^22 (S) S ) 

S-4M2 
s' - 4M2 


1 _ s' -ql- ql 

s'- s Ai2(s') 
s' s- AMl Ai2(s) / 1 

S s' — 4M2 Ai2(s') \s' — s 

_JL_ s' -ql-ql \ 

- s Ai 2 (s') ) ’ 


[7 s-4M^ Ai2(s) 1 

V s s' — 4M2 Ai 2 (s') s' — s’ 

s-AMl / 1 _ s' - A 

V S s' — 4M2 \s' — S Ai 2 (s') / 


s' - 9i - 92 \ 
Ai2(s') y’ 


(2.41) 


Because we use here a different basis, the kernel functions are slightly different from the ones given in [26]: in 
order to avoid fictitious poles in q\ — g|, one has to work with the linear combinations /i 2 ,+o i /i 2 ,o+) which 
are symmetric/antisymmetric under q^ o g|. The remaining difference to [26] is a polynomial, which can be 
reabsorbed into subtraction constants. 


2.6 Pion-pole contribution 

In order to define the contribution to 7 * 7 * —>■ 7 r“'" 7 r“ due to the exchange of a single pion, we employ a dispersive 
picture. In this subsection, we demonstrate in detail that the pure pion-pole contribution in the sense of unitarity 
is exactly the Born contribution in scalar QED (sQED), multiplied by pion vector form factors. The dispersive 
picture cleanly defines this as the correct dependence on the photon virtualities qf. We also discuss the different 
meaning of Feynman and unitarity diagrams in this context. 

To keep the presentation as transparent as possible, we omit subtleties due to isospin conventions here and 
simply present the plain sQED expressions. Compared to [26], this results in an overall sign for the amplitude of 
7 * 7 * —> 7 r+ 7 r“. For HLbL, this convention has no effect, but the choice in [26] is advantageous when transforming 
to the isospin basis and reconstructing partial waves in accordance with Watson’s theorem [42]. 

We assume that the asymptotic behavior of 7 * 7 * —>■ tt+tt” in the crossed Mandelstam variables t and u 
permits an unsubtracted fixed-s dispersion relation® for the scalar functions: 


y4-(s,t,u) 




u- TT 74^2 t' -t 


1 r ^^, Dt-uW-.s) 

TT U' -U 


(2.42) 


where pf.j ^ denote the pole residues and Df.^ „ the discontinuities along the t- and rt-channel cuts. Both are 
determined by unitarity. Consider the t-channel unitarity relation: 

Im* (^e^(2Tr76^^Hpi + P 2 - qi - 92)e^'(9i)e^"*(-92)II'^''(pi,P2, 9i)) 

1 



J dk^ (n;{fcj| 7 r {p 2 ) 7 i-<l 2 , X 2 ))*{n; {h}\TT (-pi) 7 *( 9 i, Ai)), (2.43) 


where is the symmetry factor for the intermediate state jn) and we denote the Lorentz-invariant measure by 
dk := H we single out the one-pion state from the sum over intermediate states, we find 


IrntW"" = \jdk (2^)V4)(gi - p, - k){k - p^nk + P 2 r{q^F^(qj), 


(2.44) 


where F^ is the electromagnetic form factor of the pion, defined by 

(’’■■^(fc)|jem(0)k+(p)) = {k+p7 F7{{k-pf). (2.45) 

^Whether a fixed -5 dispersion relation for 7 * 7 * tttt requires subtractions or not has no influence on the pion pole. Therefore, 
also the Mandelstam representation for the pion-box contribution to the HLbL tensor, which we discuss in this article, remains 
unaffected by possible subtractions in the sub-process. 
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Figure 2: Correspondence of the pion-pole contribution to 7 * 7 * —>■ in terms of unitarity diagrams and the Born 

contribution in terms of sQED Feynman diagrams. The dashed lines in the unitarity diagrams indicate a cut line, hence 
the internal pion is on-shell. 


After performing the trivial integral, we obtain 


iql)F^{ql) 7r6{t - {q^^ql - ql^q^, - qi^ql + q^^q'^) . (2.46) 

t- and u-channel are related by pi ^ P 2 - 

= F^iqf)F^ (g|) 7r5{u - M^) {ql^q^ - q>^ql - q^q-( + q^.q^) . (2.47) 

Note that these expressions are only gauge-invariant due to the presence of the delta function: the contraction 
of the bracket in (2.46) with q'^ or is proportional to t — 

If we project the imaginary parts of onto the scalar functions, making use of the delta function, we 
obtain 


lKA,=F:{qi)F^{qi)^5{t-Mt), 


ImiA 4 = 


2 2 

s-qf-qi 


F^{qi)F:{qi),.5{t-Mi), 


Im“Ai = F^iqf)F^{ql) ^5{u - M^), 

Im“A4 = -5- ^F^{ql)F^{ql) t^6{u - M^), 

s-qi-qi 

while the one-pion contributions to the imaginary parts of the remaining scalar functions vanish. 
The pion-pole contribution to the scalar functions is therefore 


Al = -F^{ql)F^{ql) 
Al = -F^{ql)F^{ql) 


1 


1 


t-Ml u-M^ 
2 / 1 


1 


Al = Al = Al = 0. 


' s - ql - ql \t- Ml ' u - 
We compare the pion-pole contribution with the Born contribution in sQED: 

= 



+ ie {2pl - ql;){2p’l - q^)- 


1 


- Ml 


+ 2ie‘^g>^‘', 


and read off the Born values of the scalar functions: 


^Born _ _ 


4 Born _ 

^4 — — 


1 


1 


t-Ml u-Mir 


1 


1 


s-ql-ql \t- Ml u-MlJ ' 

^Born _ ^Born _ ^Born _ q 


(2.48) 


(2.49) 


(2.50) 


(2.51) 
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Figure 3: Kinematics of the light-by-light scattering amplitude. 


We find that the pion-pole contribution corresponds exactly to the sQED Born contribution multiplied by 
electromagnetic pion form factors for the two off-shell photons.^ Note that, if we think in terms of unitarity 
diagrams, we have now considered the pure pole contribution to the scalar functions. However, in terms of 
Feynman diagrams in sQED this corresponds to a sum of two pole diagrams and the seagull diagram.® It is 
important to be aware of the different meaning of a topology in the sense of unitarity and a Feynman diagram, 
see Fig. 2. As will be shown in Sect. 5, it is exactly this distinction that makes the sQED pion loop in HLbL 
coincide with box-type unitarity diagrams representing tttt intermediate states with a pion-pole LHC, although, 
in terms of Feynman diagrams, it is composed of the sum of box, triangle, and bulb topologies. 

3 Lorentz structure of the HLbL tensor 
3.1 Definitions 

In order to study the contribution of HLbL scattering to the anomalous magnetic moment of the muon, we need 
first of all a description of the HLbL tensor. The object in question is the hadronic Green’s function of four 


electromagnetic currents, evaluated in pure QCD (i.e. with fine-structure constant a = e^/(47r) = 0): 

92,93) = -^J (3.1) 

The electromagnetic current includes only the three lightest quarks: 

Jem ■= qQYq, (3.2) 

where q = (u, d, s)^ and Q = diag(|, ——|)- 

The contraction of the HLbL tensor with polarization vectors gives the hadronic contribution to the helicity 
amplitudes for (off-shell) photon-photon scattering: 

11aiA 2.A3A4 = e^H9i)e^"(92)e^*(-93)e^’‘*(A)n''''^'^(9i,92,93)- (3.3) 

For notational convenience, we define 

54 := fc = gi-I-<72 + 93- (3.4) 

The kinematics is illustrated in Fig. 3. 

We use the following Lorentz scalars as kinematic variables — these are the usual Mandelstam variables: 

s:={qi+q2f, t := {qi + q^f, u'.= {q2 + q^f, (3.5) 

which fulfill (we will take = 0 at some later point) 

4 

s -\-t + u = ^ qf =: S. (3.6) 

i=l 

Gauge invariance requires the HLbL tensor to satisfy the Ward-Takahashi identities 

{9i ,92,93,9j}n^i.Aa(9i,92,93) = 0. (3.7) 


^Therefore, the dispersive definition of the pion pole (2.49) coincides with the gauge-invariant pole contribution of the ‘soft- 
photon amplitude’ in [43]. We thank S. Scherer for pointing this out. 

^The equivalence of the pion pole and the Born term is surprising given the fact that (2.50) contains a term with , while the 
imaginary parts (2.46) and (2.47) do not. Tracing the above steps backwards, one sees that in the t- or u-channel imaginary parts 
the coefficient of is proportional to {t — — M^) or (u — M^)S{u — M^) and hence vanishes due to the delta function. 
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3.2 Tensor decomposition 

In general, the HLbL tensor can be decomposed into 138 Lorentz structures [13, 44, 45]: 

+ E E e 

fc-1,2,4 i^l,3,4 i^l,3,4 

Z=l,2,3 Z=l,2,3 fc=l,2,4 

+ E 9‘'\^9r^Ji+ E E 

*^ 2 , 3,4 2 ^ 2 , 3 , 4 2 ^ 2 , 3 , 4 

Z=l,2,3 fc=l,2,4 j = l,3,4 

+ E E ^I'jkl 

i=2,3,4 fe=l,2,4 
i=l,3,4 (=1,2,3 

138 

=:E^r^"2,. (3.8) 


The 138 scalar functions 


{sj := {n\n2,n3,nti,n5 ,n6 ,nj„nf„n?.,ni° } 


(3.9) 


depend on six independent kinematic variables, e.g. on two Mandelstam variables s and t and the virtualities gj, 
g|, (jI, and q\. They are free of kinematic singularities but contain kinematic zeros, because they have to fulfill 
kinematic constraints required by gauge invariance. The Ward identities (3.7) impose 95 linearly independent 
relations on the scalar functions, reducing the set to 43 functions. 

As we did in Sect. 2.2 for the case of 7 * 7 * —>■ tttt, we will now construct a set of Lorentz structures and 
scalar functions, such that the scalar functions contain neither kinematic singularities nor zeros. Compared to 
7 * 7 * —> TTTT, the application of the recipe given by Bardeen, Tung [31], and Tarrach [32] is much more involved. 
Again, the recipe by Bardeen and Tung does not lead to a kinematic-free minimal basis (which would consist 
here of 43 scalar functions).® Following Tarrach, we will construct a redundant set of 54 structures, which is 
free of kinematic singularities and zeros. 

In a first step, we define the two projectors 


:= - 


U 1 / 

92 9i 
91 - 92 ' 


which have the following properties: 


9^illl = 0 , 


= n„ 


jXcr _Act 9493 

-'34 -—9 ) 

93 ■ 9i 


93lli = 0 , 

VX(7 J 

t\\ TT 
-^34 


92III = 0 , 

^(T 7-34 


9 ura = 0 , 

_ TT ^ 

^12 — -L-L/j, A<75 

7-<7 CTTT _ TT <7 

^34 ^^fii'Xcr' — 5 


(3.10) 


(3.11) 


i.e. the HLbL tensor is invariant under contraction with the projectors, but the contraction of every Lorentz 
structure produces a gauge-invariant structure. Hence, we project the tensor 


n'^'^CT = If/ 

138 

E T-/J./J. Tl/' 1/ tXX' tc' (T ri 

-^12 ^12 -^34 -^34 ^^' v ' X '( 7 '^ 


2=1 

138 43 

=:y:Lr"'H. = y;Lf; 

i=i 3=1 


X<7 


(3.12) 


Only 43 of the 138 projected structures are non-zero, i.e. all constraints imposed by gauge invariance 

are already manifestly implemented. Since the projected structures are still multiplied by the original scalar 

®We use ‘basis’ in a loose terminology: as we will discuss in Sect. 3.3, a basis in the strict mathematical sense consists of 41 
elements due to two peculiar redundancies in four space-time dimensions. 
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functions no kinematic singularities have been introduced into the scalar functions. We now have to remove 
the kinematic zeros from the scalar functions by removing the single and double poles in qi ■ and 93 • < 74 , which 
are present in the structures "We adapt the recipe of [31] (cf. Sect. 2.2): 

• remove as many (qi • q 2 , qz ■ qi) double-double poles as possible by adding to the structures linear combi¬ 
nations of other structures with coefficients containing no poles, 

• if no more double-double poles can be removed in this way, multiply the structures that still contain 
double-double poles by either qi ■ <72 or <73 • q^ (the choice is irrelevant in the end), 

• proceed in the same way with double-single, single-double poles, etc. until no poles at all are left in the 
structures. 

As already mentioned, it is again impossible to avoid introducing kinematic singularities into the scalar functions 
by applying this procedure [32]. However, the only step where kinematic singularities can be introduced is the 
multiplication of the structures by <71 • <72 or q^ ■ q^ (i.e. the division of the scalar functions by these terms). This 
means that the only possible singularities are (double or single) poles in qi ■q 2 or 93 -( 74 . The precise form of these 
poles can be easily determined: they correspond to degeneracies of the obtained basis of Lorentz structures in 
the limit gi • <72 —> 0 and/or (73 • (74 —t 0. Therefore, the 43-dimensional basis has to be extended by additional 
structures, which are found by studying the null-space of the present structures in the mentioned limits. 11 such 
structures can be found. The extended generating set of 54 structures exhibits all possible crossing symmetries 
in a manifest way. 

Explicitly, the resulting representation of the HLbL tensor reads 

54 

^ (3.13) 

i^l 

where 

^ (^q^q^ _ q^ . , 

T^X'^ = (92 91 - 9i • 929^'") (91 • 94 (91 93 - 9i ' 939^"^) + 94 9i 9i ' 93 - qUiQs ' 94 ), 
j.^^Xa ^ , ^29^'") (92 • 94 (9i 93 - 9i • 939^'") + 94 92 9i ' 93 - 9i 92 93' 94 ) , 

j^^^Xa ^ , q^gP'''^ (^q^q^ . q^ - q^q^ . (^q^q^ .q^-q^q^. ^ 

rj^^^Xa ^ ^ ^ (939492 - 9^92 93 + 9^”^ (9492 • 93 - 93 92'94) 

+ 9'"^ (92 93 • 94 - 94 92 • 93) + 9 ^'' (9392 • 94 - 92 93 • 94) ) , 

= 93 (91 • 9392 • 9494 9^'" - 92 • 9391 ' 94949^^ + 94 94 (9i 92 • 93 - 92 9i ' 93) 

+ 9i • 9493 94 92 - 92 • 9494 9391 + 9i ' 9492 • 94 - 939^“") ) 

- 94 (91 • 9492 • 93939'"'^ - 92 • 9491 • 93939^'^ + 9393 (9i 92'94 - 92 9 i ' 94) 

+ 9i • 9394 93 92 - 92 • 9393 94 91 + 9i ' 9392 • 93 ( 949^'" - 949'"'") ) 

+ 93 • 94 ( (91 94 - 9 i • 949^^) (9392 - 92 • 939 '"'") - (9294 - 92 • 949 ^'") (93 9 i - 9 i ' 939 ^'") ) ■ (3.14) 

These structures satisfy the following crossing symmetries: 

= Ci2[Tr^"] = C34K"^"] = C34[Ci2[Tr^'"]] = C24[Ci3[rr^'^]] 

= C23[Ci4[rr^"]] = C24[Cl3[C34[rr^"]]] = C23[Cl4[C34[Tr^1]], 

= Ci2[Tr^'^] = C34[TT"^"] = C34[Ci2[Tr^"]] = C24 [Ci3[T4''"^"]] 

= C23[Ci4[TT"^"]] = C24[Ci 3[C34[TT"^"]]] = ^23 [Cl4[C34 [Tr^l]], 

Xo ' ^ Acj 
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(3.15) 


= C34[Ci2[rr9"^"]], 

= Ci2[T37^n = C34[T37^"] = C34[Ci2[T37^n]. 

rpf-LuX<7 _ ^ j-^n/J.I^Ao'j 

T i/2I/A<7 ^ rr-j-if-LUX(71 ^ rrTifliyX(T'\ ^ Try^/it'A^T] ] 

49 — J — “^34[-^49 J “ *-34[<-12 49 JJ’ 

where the crossing operators Cij exchange momenta and Lorentz indices of the photons i and j, e.g.^*^ 

Ci2[f] ■= i^,qi ^ q2), Ci4[/] :=/(Ai ^ o-,gi ^- 94 ). (3.16) 

All the remaining structures are just crossed versions of the above seven structures, as shown in App. B.l. 
Since the HLbL tensor is totally crossing symmetric, the scalar functions Hi have to fulfill exactly 

the same crossing properties of the corresponding Lorentz structures as given in (3.15) (note the antisymmetric 
crossing relations in 1149 ). Therefore, only seven different scalar functions Hi appear, together with their crossed 
versions. These scalar functions are free of kinematic singularities and zeros and hence fulfill a Mandelstam 
representation.^^ They are suitable quantities for a dispersive description. 

The subset consisting of the following 43 Lorentz structures forms a basis: 

e {1,...,21,23,25,27,29,30,33,...,36,38,...,45,49,...,53}I. (3.17) 

The corresponding scalar coefficient functions 11^, defined by 

43 

^ (3.18) 

i=l 

exhibit kinematic singularities in qi ■ and qs ■ q^. The exact form of these kinematic singularities can be 
determined by projecting (3.13) on this basis. The relation between the basis coefficient functions 11^ and the 
scalar functions Lti is given in App. B.2. 

While crossing symmetry is manifest if the HLbL tensor is expressed in terms of the redundant set of 54 
structures, it is partially obscured in the basis of 43 elements. Still, many basis functions 11^ are related by 
crossing symmetry, see (B.3) in App. B.2. The basis coefficient functions are completely specified by crossing 
symmetry and 9 representative functions, e.g. 

|f[l, 144 , fty, ft9, ni9, ft21, Hag, f[39, n4o| . (3.19) 

We stress again that our decomposition (3.13) is both manifestly crossing symmetric and gauge invariant. In 
particular, this implies that all soft-photon zeros [39] of the helicity amplitudes (3.3) are already implemented, 
in contrast to the helicity basis employed in [26]. In this case, only specific combinations of basis functions 
were free of kinematic singularities and zeros and thus admitted a dispersive representation, while the soft- 
photon constraints needed to be imposed by hand. The BTT set (3.14) provides a Lorentz decomposition 
that incorporates soft-photon constraints and avoids kinematics by construction. On the other hand, the main 
motivation for choosing the helicity basis was that the unitarity relations in the form of helicity partial-wave 
amplitudes become diagonal. Of course, this will not be true in the BTT basis, so that the derivation of an 
explicit dispersive representation for the BTT amplitudes requires a basis change. In practice, this change 
between BTT and helicity bases proceeds by means of (3.8) as an intermediate step: the helicity amplitudes 
can be easily expressed in terms of (3.8), while the remaining basis change can be performed along the lines 
explained in App. C. 

3.3 Redundancies in the Lorentz decomposition 

Because the set of 54 scalar functions is redundant, there is obviously an ambiguity in the definition of the coef¬ 
ficient functions. However, crossing symmetry and the absence of kinematic singularities restrict this ambiguity 
in a very specific way. The shifts 

Hj l-^> Hi -I- JHi (3.20) 

^®The composition of two crossing operators is understood to act e.g. in the following way: Ci 2 [C 23 [/(<?i, <i' 2 , <? 3 , 94 )]] = 
Ci2[/(9l,g3,<?2,94)] = /(92,93, 91, 94 ). 

^^In principle, the crossing antisymmetry of 1149 implies a kinematic zero. We ignore it here and show in App. F.3 that this zero 
has no impact on the dispersion relation for the HLbL tensor. 
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leaving the HLbL tensor invariant are easily found by requiring that the basis elements 11^ remain 

invariant. We find 

6Ui = 0 , sn^ = 0 , 

(5117 = —q2 ■ q3q2 ■ (Z 4 A 19 , ^Ilig = qi ■ < 74(72 • <73^19: ^Ilai = —q^ ■ g4Ai9, (3.21) 

51137 = <71 • (72A37, 51149 = <71 • (72A49, 


where the shifts have to satisfy 


^19 = Ci 2 [Ai 9 ] = C 34 [Ai 9 ] = C 34 [C 12 [A 19 ]], 

A 37 = C24[A37] = C23[A37] = C34[A37] = C24[C34[A37]] = C 24 [C 23 [A 37 ]], 

A 49 = —Ci 2 [A 49 ] = —Ci 3 [A 49 ] = —Ci 4 [A 49 ] = —C23[A49] = —C24[A49] = —C 34 [A 49 ] = . . . . (3.22) 


The ellipses in the last line stand for the permutations that are not listed explicitly: A 49 is symmetric (anti¬ 
symmetric) under all even (odd) crossing permutations. The shifts 511^ that are not specified follow from the 
crossing relations (B.l). 

Apart from the redundancy introduced by extending the basis to a set of 54 elements, there exists another 
ambiguity in the definition of the scalar function. As it was pointed out in [46], the 138 initial Lorentz structures 
of the HLbL tensor fulfill two linear relations in d = 4 space-time dimensions. This is not obvious in a covariant 
notation, but can be easily shown in an explicit reference frame. The two linear relations can be translated into 
two linear relations between the 43 basis elements (which means that these structures actually do not form a 
basis in a strict sense). Therefore, in four space-time dimensions there are only 41 independent scalar functions, 
which matches the number of fully off-shell helicity amplitudes. This is analogous to the case of 7 * 7 * —t tttt, 
where the number of helicity amplitudes (5) also matches the number of independent scalar functions. 

The two linear relations can be further translated into relations amongst the 54 elements of the redundant 
crossing-symmetric set of Lorentz structures. In fact, it turns out that expressed in this way, the two relations 
become particularly simple. They are given by: 

54 

= 0, (3.23) 

i=l 

where 


c\ = —s -\-1 + u, c\ = —s — t + u, C 4 = —5s — t + u, C 5 = 2(—s -|- m), 

„i _ „i _ „i _ „i _ „i _ „i _ „i _ „i _ „i _ „i _ 4 

^7 — ^8 — '-17 — ''18 “ ''19 ~ ''20 “ ''29 ~ ''30 “ ''50 ~ ''53 “ 

pl - A - pi - A - - pi - _4 

''ll — '^12 — ''13 ~ ''14 — ''23 “ ''24 — ''25 ~ ''26 ~ 

''37 = ''46 = ''47 = ''48 = ''51 = ''52 = 2, 

''38 = ''40 = ''43 = ''44 = '^49 = ''54 = 


= —s +1 + 5u, 


(3.24) 


and 


C2 = S — t-l-M, C3 = —S — t-|-M, C4 = 2 ( —t-I-u), C5 = —s — 5 t 4 

22222222.1 
'^9 ~ '^10 ~ '^15 ~ '^16 “ '^21 — '^22 — '^27 ~ '^28 ~ 4, 

2 _ 2 _ 2 _ 2 _ 2 _ 2 _ 2 _ 2 _ 2 _ 2 _ 2 _ 4 

'^11 — '^12 — '^13 ~ '^14 — '^14 — '^23 ~ '^24 — '^25 “ '^26 ~ '^49 ~ '^54 — 4 


''38 


''39 


_ 2_2 _ 2 _ 2 _ 2 _ n 

— '^40 ~ '^43 ~ '^44 — '^51 ~ '^52 — 

2 2 2 2 2 0 

— — C42 — 045 — Crg — C53 — 2 , 


= s — t + 5u, 


(3.25) 


and where all unspecified cj’^ are zero. Eq. (3.23) is most easily verified by contracting the equation with the 
138 tensor structures. We have checked explicitly that the matrix is of rank 136. The linear 

relations between the structures imply the following ambiguities in the scalar functions: 

n, ^n, + ci/i + c2/2, (3.26) 


where /i ^2 are a priori unspecified functions free of kinematic singularities. The requirement that the scalar 
functions satisfy crossing relations of the same form as (3.15) and (B.l) imposes crossing relations on the 
functions fi^ 2 - 

fl = Cl 2 [/ 2 ] = C34[Ci2[/l]] = C24[Ci3[/i]] = Culfi] = —Ci 2 [/l] — Ci3[/i]. (3.27) 
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Let us write 


1149 = ^ (1149 + C24[Cl3[n49]] +Ci4[n49] + C14[C24[Cl3 [1149]]]) 

+ ^ (1149 + C24[Ci3[n4g]] — C14P49] — Ci4[C24[Ci3[n49]]]) 

+ 2 (^49 — C24[Ci3[n49]]) . (3.28) 

The first bracket in this expression is given by 1149 — IIso — 1153 +1154 and has exactly the crossing properties 
of fi- Therefore, the ambiguity can be fixed by the condition 


1149 — IIso — 1153 + 1154 — 0. (3.29) 

Very recently, a Lorentz decomposition of the HLbL tensor into 72 redundant structures has been presented 
in [47]. A comparison with our redundant set of 54 structures, as first given in [33], has been performed, 
confirming the absence of kinematic singularities. By studying systematically the irreducible representations 
of the permutation group 5 * 4 , the authors of [47] aim at constructing a minimal basis consisting of 41 totally 
crossing symmetric elements free of kinematic singularities. Note, however, that [47] neither provides this basis 
nor a concrete recipe for how to obtain it. Therefore, it remains to be seen if a minimal basis free of kinematic 
singularities exists at all. For our purpose the question whether such a basis exists is immaterial, since for a 
dispersive representation of the HLbL tensor all that is required is a decomposition into scalar amplitudes free 
of kinematic zeros and singularities. 

To conclude this subsection, we stress that the presence of kinematic singularities in any known basis requires 
us to work with a redundant set of 54 Lorentz structures and scalar functions H^. Although the redundancy 
of the Lorentz structures introduces an ambiguity in the functions H^ (which is already very restricted by 
crossing symmetry and will be further restricted by analyticity), this ambiguity will cancel in the calculation 
of any physical quantity such as {g — 2)^. In particular, this implies that the additional redundancy in d = 4 
dimensions does not affect our calculation in any way, since the minimal set suited for a dispersive representation 
still involves 54 functions. In the rest of the paper we will continue to work with the BTT amplitudes defined 
in (3.14). 


4 HLbL contribution to {g — 2)^ 

In Sect. 4.1, we review the definition and calculation of a^, as well as general techniques for the calculation of 
(see e.g. [48]). The well-known general formula requires still a rather long calculation before a number 
can be finally obtained. For the pion-pole contribution, these steps have been worked out long ago [19]. With 
our complete set of 54 kinematic-free structures (3.13), this procedure can be repeated for the whole HLbL 
contribution in full generality, as we explain in Sects. 4.2 and 4.3. In Sect. 4.4, we study the validity of the 
applied Wick rotation in the presence of anomalous thresholds. 

4.1 Projector techniques 

The T-matrix element of the interaction of a muon with the electromagnetic field is defined by 

T''(pi,P 2 ) := (4.1) 
where e = jej. Diagrammatically 



Assuming parity conservation, the vertex function F^ can be decomposed into form factors as 

r^{pi,P2) = rF,{e) - ^^^F2{k^), 

2 ? 77 .^ 


(4.3) 
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where k = pi — P 2 and := §[ 7 ^, 7 '^]- Fi is the electric charge or Dirac form factor, F 2 the magnetic or Pauli 
form factor. The anomalous magnetic moment of the muon is given by 

a^ = \{g-2)^ = F2{Q). (4.4) 

We expand the vertex function T'' to first order in powers of k^: 


r''(pi,P2) = ^^{p,p) + k^—T^(pi,p2) 
=: P''(p) + A:,r^"(p) + ..., 




(4.5) 


where p := \{pi + p 2 )- 

Using projector techniques together with angular averages (see [19, 48]), the anomalous magnetic moment 
of the muon can be calculated as 


= Tr 


48mu 


1 

12 ' 

Tr( 




m,, 


1 p^^ 

4 ma 


' + w^)[7^, 7l(|> + mf,)rf,p{p)) 


where now p^ = 

We are interested in the contribution of the HLbL tensor to a^, diagrammatically 


where 



1 


9i92(pi -P2-qi 


(P2)r()LbL(Pl.P2)M(pi), 


(j ^2 +W^) (^1 - ^2 

(P2 + gi)2 - to2 (pi - q2Y - rrijP'' 

{qi,q2,pi -p2-qi- 92 ). 


(4.6) 


(4.7) 


(4.8) 


The HLbL tensor is defined in (3.1). Differentiating the fourth Ward identity in (3.7) with respect to kp = 
{qi + q 2 + q 3 )p yields 

d 

n^i.Ap(gi,92,fc - <?i - 92 ) = -k'^-^^p,,\a{qi,q2,k- qi - 92 )- (4.9) 


It was already argued in [49] that vanishes linearly with k (i.e. the derivative contains no singularity), and 

so must This is easily verified with our tensor decomposition (3.13). Therefore, the HLbL contribution 

to the anomalous magnetic moment is given by 


,HLbL 


48m, 


-Tr ((|^ + mp)[Y, 7‘"](|> + , 


(4.10) 


where 




fc =0 


We use the Ward identity (4.9) to write 

rHLbL.„ ,6 [ d% dU2 {p2 + .7 - ^2 + 

^ ’ J ( 27 r )4 ( 27 r) 4 ^ (p 2 +gi) 2 -my (pi - 52 )^ - 


1 


qfqUpi -P 2 -qi- 92)^ dkp 


d 

k'^^r-IlpuXaiqi,q2,k- qi - < 72 )- 


^Note that we have defined k as outgoing, resulting in the different sign of the second term with respect to [48]. 


(4.11) 


(4.12) 
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Taking the derivative and limit leads to the well-known expression 


pHLbL 

■*- p <7 


J ( 27 r )4 ( 27 r) 4 ^ (p + g^)2_^2 ' 
1 9 

qMiqi + <72)2 52, fc - 


A 


+w^) 
(jp-q2Y - ml 


qi - 92) 




(4.13) 


4.2 Loop integration 

In order to compute the contribution to {g — 2)^, one has to take the trace in (4.10) and perform the two-loop 
integral of equation (4.13). Five of the eight integrals can be carried out analytically with the help of Gegenbauer 
polynomial techniques [50]. To this end, we employ the representation of the HLbL tensor in terms of the 54 
Lorentz structures 


HLbL 


e® f d'^qi d'^q 2 111 

48m^ J (27r)4 (27r)4 qfq^iqi + 92 )^ (p + gi)^ - ml {p - q2)‘^ - ml 

X Tr + mi,)[Y + + 


54 

E 

i=l 


-^Ti^xaiqi,q2,k- qi -g2) 


n*(gi,g2,-gi - g2)- 


k^O 


(4.14) 


It turns out that there are only 19 independent linear combinations of the structures that contribute to 

{g — 2)^. It is possible to make a basis change in the 54 structures 

54 54 

^ Tf(4.15) 

i=l i=l 


in such a way that in the limit k —>■ 0 the derivative of 35 structures vanishes. Although this change 

of basis does not introduce kinematic singularities into the scalar functions Hi, it somewhat obscures crossing 
symmetry. The 19 structures that contribute to {g — 2)^ can be chosen as follows: 


J = ^ igj = 


i= 1,...,11,13,14,16,17 


} 


U rp^uXcT rp^uXa rjifll^Xa rp^luXcT ^^0^ 


The corresponding 19 scalar functions Hi are linear combinations of 33 scalar functions Hi, see App. D. The 21 
scalar functions not involved in these relations 


|ni i = 12,15,18,23,24,27,28,29,30,32,35,36,37,38,41,44,45,48,49,52,53| 


are irrelevant for the calculation of {g — 2 )^. 

The HLbL contribution to {g — 2)^ can now be written as 


where 


„HLbL 


—e 


d'^qi d^q2 111 

(27r)4 (27r)4 qlql{qi + ga)^ {p + gi)^ - ml {p - ga)^ - 
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x5:T,(gi ,q2;p)fit{qi,q2,-qi - g2), 

i=l 


T^{ql,q2;p) ■■= ^^^Tr ((^-f m^)[7'’,7'"](|i-kTO^)7''(^-f + to^)7 + mf,)j'- 


d - 

X ( ^T))i.A<T(gi,g2,fc-gi - g2) 


fc =0 


(4.17) 


(4.18) 


(4.19) 
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and the 11^ are needed for the reduced kinematics 

s={qi + q 2 f, t = ql, u = ql, qj, ql, ^ (gi + g 2 )^ k'^=qj=0. (4.20) 

The explicit result of the trace calculation and the contraction of the Lorentz indices is given in App. E.l. 

We can reduce the number of terms contributing to (g — 2)^ further by using the symmetry under the 
exchange of the momenta qi o —q 2 '. the loop integration measure and the product of propagators are invariant 
under this transformation, while the kernels transform under qi ^ —q 2 as 

rji . . rji rji . . rri rri . . rri rri . . rri 

Jl -r ±2 ' ' -^3, J4 -K - r J4, J5 i - r Jg, 

rji . . rri rri . . rj~\ rT~i . . rj~\ rj~\ . . rj~\ 

It i - > J-8^ ^ - > -L12: ^ 10 ^^11^- ^^14) 

Ti5 <—)• Ti5, Ti 6 <—>■ Tie, Tit ^^ ^ig, Tig <—>■ —Tig. (4.21) 


For the reduced kinematics (4.20) the exchange qi O — (72 is equivalent to the crossing transformation t ^ u, 
qf o q 2 . With the help of the crossing relations of the scalar functions 11^, it is easy to check that the Hi 
transform analogously to the kernels Tj, i.e. 

Hi ^—> Hi, fig ^—> Hg, n4 ^—> 114, rig <—)■ fig, 

fly ^- > Hg, fig < > fil 2 , flio ^^ riig, Hu ^- 5 - fli4, 

His ^^ hlis, rile <—)■ ni6 ) Hit ^^ hlig, Hig <——flig. (4.22) 

Therefore, it is convenient to write the HLbL contribution to {g — 2)^ as a sum of 12 terms: 


_ 6 f <^'‘gi d‘^q2 _1_1_1_ 

® J (27r)4 (27r)4 qjq'^iqi + < 72 )^ (p + gi)^ - {p - 92 )^ - 
12 

X X!^i^»i(9i,g2;p)n»,(gi,g2,-gi - 92 ), (4.23) 


where 


{ij I j = 1,..., 12} = {1,2,4,5,7,9,10,14,15,16,17,19}, 

le, |j = 1,..., 12} = {1, 2,1, 2,2, 2, 2, 2,1,1, 2,1}. (4.24) 

Note that the first two terms in this sum correspond to the well-known result for the pion-pole contribution [19] 
(up to conventions: exchange of Ti and T 2 , the explicit factor ^2 = 2, and symmetrization of Ti). 

In (4.23), the integrand depends on the five scalar products q\, g|, qi-q 2 , p-qi, a,ndp-q 2 , where the dependence 
on the last two is given explicitly (the scalar functions only depend on q^, g|, and qi ■ ( 72 )- Therefore, five of 
the eight integrals can be performed without knowledge of the scalar functions. The same integrals as in the 
case of the pion-pole contribution occur [3, 19], which have been solved with the technique of Gegenbauer 
polynomials ]50]. This method has been applied before to the full HLbL contribution in the context of vector- 
meson-dominance and hidden-local-symmetry models [51, 52]. 

We perform a Wick rotation of the momenta < 71 , q 2 , and p (see Sect. 4.4) and denote the Wick-rotated 
Euclidean momenta by capital letters Qi, Q 2 , and P. Note that Q\ = —qi, Q\ = —(zli Since 

is a pure number, it does not depend on the direction of the momentum P of the muon, hence we can take the 
angular average by integrating over the four-dimensional hypersphere: 

HLbL _ 

— 


d^4:{P) HLbL 
2^2 ■ 


(4.25) 
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The kernels are at most quadratic in p, therefore we need the following angular integrals [3]: 
f (i 04 (P) 1 11 


1 1 

Y +ml 

{P-Q2Y+ml 

dn^p) 

1 

2tt'^ 

{P + Qi^ + ml 

dn^p) 

1 

2tt'^ 

{P-Q2Y+ml 

dn^p) 

P-Q2 

2tt'^ 

{P + QiY + ml 

dn^p) 

P-Qi 


-arctan 


zx 

1 — ZT 


2^2 (p - + ml 


1 -^f 

2 m 2 

1 -^ 2 ^ 

2 m 2 


— Qi' Q 2 


{l-aff 




( 1 ^)^ 

8 to ?, 


— —Qi • Q2 


(4.26) 


where r = cos 6 * 4 , defined by Qi-Q 2 = |( 5 i||Q 2 |t, is the cosine of the angle between the Euclidean four-momenta 
Qi and Q 2 , and further 


/ 4m^ 

o-f:=Wl + -^, ^12 := |( 3 i||Q 2 k, a; := Vl - r 2 

_ IQ1IIQ2I „E\ 

^ .2 ^2 )■ 

4to^ 


(4.27) 


4.3 Master formula 

After using the angular integrals (4.26), we can immediately perform five of the eight loop integrals by changing 
to spherical coordinates in four dimensions. This leads us to a master formula for the HLbL contribution to the 
anomalous magnetic moment of the muon: 


HLbL 


2o? 

Stt^ 


^00 ^00 pi _ 12 

/ dQi / dQ2 / dT\/T^^Q\Ql'^Ti{Qi,Q2,T)Iii{Qi,Q2,T) 
Jo Jo J-i ,■ 1 


(4.28) 


where Qi := |(5i|, Q 2 ■= \Q 2 \- The hadronic scalar functions 11^ are just a subset of the Hi and defined in (D.2). 
They have to be evaluated for the reduced kinematics 

s =-Qi =-Qi - 2 Q 1 Q 2 T - Ql, t = -Ql, u = -Ql, 

= ql = -Ql, ql =-Ql =-Ql-‘J‘QiQ 2 T-Ql, = o. (4.29) 

The integral kernels listed in App. E.2 are fully general for any light-by-light process, while the scalar 
functions 11^ parametrize the hadronic content of the master formula. In particular, (4.28) can be considered 
a generalization of the three-dimensional integral formula for the pion-pole contribution [3]. It is valid for the 
whole HLbL contribution and completely generic, i.e. it can be used to compute the HLbL contribution to 
{g — 2)f^ for any representation of the HLbL tensor, irrespective of whether the scalar functions are subsequently 
specified dispersively or taken from a model calculation. For an arbitrary representation of the HLbL tensor, 
the scalar functions H^ can be obtained by projection, see App. C and App. F.2. 

An analogous master formula was derived in [26] in the case of a helicity basis for the HLbL tensor. This step 
in the calculation is completely equivalent, in both cases the calculation of proceeds via an evaluation 

of the trace in (4.10) and subsequent reduction of the two-loop integral with Gegenbauer techniques. From a 
technical point of view, the BTT approach offers several simplifications, since the H-wave-related amplitudes 
in [26] require another angular average to define the A: —)■ 0 limit as well as more complicated Gegenbauer 
integrals than the standard ones given in (4.26). 

In close analogy to the pion-pole contribution [19], the main benefit of the master formula (4.28) is the fact 
that it contains only a three-dimensional integral, and thus is well-suited for a direct numerical implementation. 
In particular, the energy regions generating the bulk of the contribution can be identified by numerically 
integrating over t and plotting the integrand as a function of Qi and Q 2 [19, 51-53]. 

Before turning to the main part of this paper, the foundations for a model-independent calculation of the 
scalar functions H^ by making use of dispersion relations, we next consider a subtlety in the application of the 
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Gegenbauer integrals in (4.26). To obtain these integrals a Wick rotation has been performed, whose validity 
might be questioned if the HLbL tensor possesses non-trivial analytic properties. For instance, already a simple 
triangle loop function would give rise to anomalous thresholds due to the integration over all photon virtualities, 
see [54-56]. In the following subsection, we prove that the Wick rotation remains valid even in the presence of 
anomalous thresholds. 


4.4 Wick rotation and anomalous thresholds 


Performing a Wick rotation before calculating a loop integral is a standard procedure. The applicability of this 
transformation depends on whether the integrand is free of singularities in the region swept by the integration 
path during the rotation. In our case we have to perform the two-loop integral (4.23), with an integrand 
containing products of propagators multiplying functions satisfying a Mandelstam representation, so having 
the analytic properties of one-loop functions. The case of the propagators is standard. The one-loop functions 
occurring in our integrand, however, need some more discussion. The complication is not due to the presence of 
cuts instead of poles: cuts of loop functions are usually on the positive real axis above some threshold, and the 
functions are evaluated on the upper rim of the cut. The Wick rotation can be applied without problem even 
in these cases. The need for a more detailed discussion is due to the presence of anomalous cuts — cuts which 
appear first in three-point functions and which, depending on the value of the momenta squared of the external 
legs, may intrude into the first Riemann sheet away from the real axis. In the present paper we need to deal 
with three- and four-point one-loop functions, such as C'o(s;gj,g|) and Do(s, t; ( 7 i, gf, 94 ) (see App. G). For 
a discussion of the analytic properties of these loop functions, see [55, 57-60]. We will illustrate what happens 
in the presence of such anomalous cuts by discussing the case of the Co function: the conclusion is that the 
Wick rotation can be applied without problem even in these cases. 

The function Co{s\q\,q2) has an anomalous cut for q\ > Q, q^ > 0, and q\ + q2 > 4to^, where m is the 
mass of the particle running in the loop (assuming it is only one mass). The origin of this anomalous cut is 
that Co admits a dispersive representation having a discontinuity on the positive real axis for s > 4m^. The 
discontinuity, however, has itself a cut between the branch points: 

s± = 9? + <zl - T - 4m2)(j|((72 - 4m2). (4.30) 


While for qf 2 < 0 these branch points lie on the second Riemann sheet, for qf > 0, > 0, and qf+q2 > 4to^ one 

of them moves into the first Riemann sheet by crossing the real axis above s = Am?. In this case the dispersive 
representation of the function Cq needs to be modified by adding a further integral along the anomalous cut 
from 4 to^ up to s+. The most convenient way to see how the branch points move as a function of qf and q^ is 
to introduce the variables [55] 




(4.31) 


The upper half-plane of the complex variables qf are mapped onto the upper unit semicircle of unit radius: the 
negative real axis of qf is mapped onto the segment 0 < < 1 , the region 0 < qf < 4m^ onto the semicircle 

with 0 < (j)i < TT, and finally the semiaxis 4m^ < qf onto the segment — 1 < < 0. The branch points 

s± can be expressed as follows in terms of these variables 


_ (1 - a )^ (1 - 6 )^ (1 - gi)^(i - 6 )^ (1 - ei^)(i - g|) . , 

Cl C 2 2 C 1 C 2 2 C 1 C 2 

and from this it is easy to check that s_ always stays on the second sheet, whereas s+ moves into the first one 
under the conditions given above. It is important to stress that when it enters the first Riemann sheet, s+ stays 
always below the real axis, either on its lower rim (for qf < Am? and qf < 4m^ or for qf > 4m^ and qf > 4to^), 
or into the lower complex half-plane (for qf < Am^ and qf > 4to^ or vice versa), in contrast to what is shown 
in Fig. 2 of [55]. 

When we perform the Wick rotation for qi and <72 (to avoid weird trajectories in the complex s-plane we 
rotate both simultaneously), their squares (in Minkowski space) become transformed into minus their Euclidean 
squares {—Qf), and can only take negative values. But also s = (gi -I- < 72 )^ is mapped onto se = —{Qi + ( 52 )^, 
which equally can only be negative. If we start from a positive value of s and perform the Wick rotation, 
the trajectory of s follows an arc in the upper complex half plane and lands on the negative axis at Se- The 
potential danger is if we start from s > 4m^, just on the upper rim of the cut, but the Wick rotation moves it 
away both from the normal as well as from the anomalous cut, which, as mentioned above, is always below the 
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Figure 4 : Movement of the trajectory of the anomalous threshold s+ during the Wick rotation. The thick gray line on 
the real axis is the unitarity cut for s > 4 m^. The solid line denotes the trajectory of for a fixed qf > 0 : the anomalous 
threshold enters the physical sheet through the unitarity cut when 52 increases. The dotted lines show the movement 
of when qi and q2 are simultaneously Wick rotated: it retreats through the unitarity cut into the second Riemann 
sheet. The dashed lines show intermediate trajectories of s+ during the Wick rotation. The part of the trajectories that 
lies on the second sheet is shown in light gray, the part on the physical sheet is black. 


real axis. If we start from negative values of s we also follow an arc in the upper complex half plane and land 
at se < s, without ever coming even close to the cuts. 

The discussion above is still incomplete, however, because we have so far neglected the fact that as we make 
the Wick rotation not only s moves, but also the qf and with them s_|_, which determines where the anomalous 
cut is. In principle, S-|- could move into the complex upper half-plane and interfere with the trajectory of s 
during the Wick rotation. An explicit calculation shows, however, that this does not happen: at the end of the 
Wick rotation both qf have become —Qi and therefore negative, which implies that the positivity conditions 
required for the existence of the anomalous cut are not fulfilled anymore. The trajectory of s_|_ during the Wick 
rotation is shown in Fig. 4. There one sees that in all possible cases s+ moves through the unitarity cut back 
into the second sheet, so that indeed, the anomalous cut disappears. The two figures show the two possible 
cases, where the anomalous threshold enters the physical sheet. Fig. 4a shows the trajectory of the anomalous 
threshold for a fixed value of 0 < < 4m^ and g| varying from —00 to 00 (solid line, light gray indicates the 

part of the trajectory on the second sheet, whereas the part on the physical sheet is shown in black). Fig. 4b 
shows the analogous trajectory for a fixed value of q\ > 4m^. The dashed lines are intermediate trajectories, 
when the two momenta are simultaneously Wick rotated by the same phase. The dotted lines indicate the 
movement of selected points on the trajectory of the anomalous threshold during the Wick rotation. All points 
of the trajectory that originally lie on the physical sheet move back through the unitarity cut into the second 
sheet. 

As a further check we have calculated numerically a simplified toy example that exhibits all relevant features 
discussed above. We consider 




d^qi 

{2eY 


__ 

(27r)4 ((^2 _|_ iYj (g| -I- ie) (s -I- ie) (s - -|- ie) 


dx 

lo Jo 


dyA{x,y) 


-1 


(4.33) 


where 

A(a;, y)= Ml- xys - x{l - x - y)ql - y{l - x - y)ql - ie. (4.34) 

The function defined by the last two integrals is proportional to the Feynman-parameter representation of the 
triangle loop function C'o(s;gi,g|), while the remainder of is chosen in as close analogy to the real HLbL 
integral as possible, including a sufficient number of propagators to make the integral converge. In the standard 
loop-integral approach the intricacies due to anomalous thresholds are automatically taken into account by the 
le-prescription in (4.33) and (4.34), leading to the following Feynman-parameter representation 


= - dx 


dy / dxi 


dX2 


dX3 




Ai{xi,y)A2{xi,X2, X 3 , X, yY 


Ai{xi,y) = l-xi+ xiy{l - y), 
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(4.35) 


+ xixy)"^ 

iixi,y) 

In analogy to the HLbL master formula (4.28) the Wick-rotated representation becomes (the Qi have been 
rescaled by m^) 


A2{xi,X2,X3,x,y) = (l - X3)Ai{xi,y) -I- X3(x2 + a;ia;(l - x)) - X3 


(X2 


A 


ar = -- rdQ^ r dQ2 f f dx f \yA{x,y)-\ 

Jo Jo J-1 Qi{^ + Qi) Jo Jo 

J^ix, y) = ^ + xyQl + x{l-x- y)Q‘l + y{l-x- y)Ql. (4.36) 

ml 

We checked numerically that indeed the representations (4.35) and (4.36) are identical. This confirms that the 
Wick rotation is permitted despite the occurrence of anomalous thresholds, all signs of which disappear in the 
Feynman parameterization as long as the le-prescription is applied consistently. Similarly, the final result for 
the Wick rotation displays no remnants of the anomalous thresholds. As shown in Fig. 4, in this case the deeper 
reason can be traced back to the trajectory of the anomalous threshold during the Wick rotation towards the 
second sheet. 

Finally, we note that the above discussion of the anomalous threshold of the triangle diagram is already 
sufficient for the full HLbL calculation. Box diagrams do not lead to further complications, since, in the master 
formula (4.23), the limit k —>■ 0 is taken before the Wick rotation of the loop momenta, hence the analytic 
structure of box diagrams is already reduced to the one of triangle diagrams. 


5 Mandelstam representation 

In the previous section, we have obtained a master formula (4.28) for the HLbL contribution to the anomalous 
magnetic moment of the muon, where the hadronic dynamics is parametrized in terms of the scalar functions 
Hi. Since these functions are free of kinematic singularities and zeros, they are the quantities that should satisfy 
a Mandelstam representation [61]. We need to determine seven scalar functions that are not related to each 
other by crossing symmetry. Due to the complexity of the problem, we cannot obtain an exact solution for the 
scalar functions but have to rely on approximations. Our strategy is to order the contributions to the scalar 
functions by the mass of the intermediate states. The lightest intermediate states are expected to give the 
most important contribution, heavier states are suppressed by the higher threshold and smaller phase space. 
In this paper, we will consider the two lowest-lying contributions: one- and two-pion intermediate states in all 
channels, i.e. pion-pole and pion-box topologies. These two contributions can be described without any further 
approximation. Of course, pion-pole and pion-loop contributions have been considered before in various model 
calculations. The important point of the dispersive approach is that these contributions become unambiguously 
defined: the pion pole, with on-shell pion transition form factor, corresponds precisely to one-pion intermediate 
states, and the sQED pion loop dressed with pion vector form factors to two-pion intermediate states with pion- 
pole LHC. The explicit proof of these identifications will be presented in this section, based on a Mandelstam 
representation for the BTT scalar functions. 


5.1 Derivation of the double-spectral representation 

For the derivation of a Mandelstam representation of the scalar functions, we follow the discussion in [62]. We 
assume that the photon virtualities qf are fixed and small enough so that no anomalous thresholds are present. 
A parameter-free description of the HLbL tensor and therefore an unsubtracted dispersion relation is crucial. 
The behavior of the imaginary parts, which is determined by the asymptotics of the sub-processes, does indeed 
suggest that no subtractions are needed. Furthermore, even a quark-loop contribution to the HLbL tensor has 
an asymptotic behavior that requires no subtractions.^^ Hence, for a generic scalar function H^, we write a 
fixed-t dispersion relation without any subtractions: 


n‘(s,t,u) =4 + 


Pi:s 


s- Ml 


Pi;u 

U — Ml 


^ 1 /■“ 

JiMl s' — S 


^ JiMl U' — U 


(5.1) 


^^Contrary to possible subtractions in the sub-processes, the presence of subtraction constants in the HLbL scalar functions 
would imply a contribution to that is not determined by unitarity. 
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Figure 5: Intermediate states in the direct channel: pion pole and two-pion cut. 


where c* is supposed to behave as lime* = 0 and takes into account the t-channel pole. The imaginary parts 
are understood to be evaluated just above the corresponding cut. The primed variables fulfill 

4 

s' + t + u' = Y.-=^qj. (5.2) 

i=l 

If we continue the fixed-t dispersion relation analytically in t, we have to replace the imaginary parts by the 
discontinuities, defined by 


hence 


Di-sis') ■= ^ i^lis' + ie,t,u') - n‘(s' - ie,t,u')) , 
Dl;uiu') ■= Yi + *e) - Illis',t,u' - ie)) , 


nl{s,t,u) = 4 + 


P%s 


s-Ml 


+ 


--Ml 


r°° D* 

/ ds'-^ 


S' — S TT . 


du 


,DUW) 




(5.3) 


(5.4) 


Both the discontinuities as well as the pole residues are determined by s- or ■u-channel unitarity, which also 
defines their analytic continuation in t. While „ are due to a one-pion intermediate state, „ are due to 
multi-particle intermediate states, see Fig. 5. We limit ourselves to two-pion intermediate states and neglect 
the contribution of heavier intermediate states to the discontinuities. 

First, we study the pion-pole contribution by analyzing the unitarity relation: 


Im^ (e*{2Tr)'^S^'^4qi + q2 + qs - qA)HxiX2MM^ 

= in /(7*(-93,A3)7*(<?4,A4)|n;{pJ)*(7*((7i,Ai)7*(g2,A2)|n;{pi}), (5.5) 


where 5'„ is the symmetry factor of the intermediate state \n). We consider now only the 7r° intermediate state 
in the sum: 

ImJ (^e'^{2TT)'^SD\q^ + ^2 + ga - q4)HxiX2,\3>^i) 

= \ J (7*(-g3,A3)7*(g4,A4)|7r°(p))*(7*(gi, Ai)7*(g2,A2)|7r°(p)). (5.6) 

After reducing the matrix elements and using the definition of the pion transition form factor 

i y e*« "^(0|r{j^^(a;)je^,„(0)}|7r°(p)) = (g2, (g - p)^), (5.7) 


we find 




-I I dp (27r)4d(4)(gi +q,- p)e*‘^“^e^'^^^gi„g2^g3W45 

X (^qy , q2) (g3: g4) 

-7rd(s - Ml)e'"'°'^qi^q2pq^.yqAsd^7T0^>^* (g?, {ql,ql) ■ 


(5.8) 
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By projecting onto the scalar functions Hj, this leads to 

t / *^ 7707 * 7 * (^1 j 92)‘^7r07*7* (^3 7 ^ 4 ) 

Pi-S - ^ Q 

and, analogously, 

t _ / *^ 7707 * 7 * (<?i , ^?4)*^7r07*7* (^2 7 ^ 3 ) 

ri\u I Q 

In order to identify the discontinuities, we project the unitarity relation selecting two-pion intermediate 
states: 

Inis’' + (?2 + <73 - 

= \ J ^=^Pi^^P2(7i’^(Pi)7i”(P2)|7*(-93, A3)7*((74, A4))*(7r+(pi)7r"(p2)|7*(<7i. Ai)7*(92, A 2 )) 

+ \I '^PWP2(7r°(Pi)7i°(7i2)|7*(-<73, A3)7*(94, A4))*(7r°(pi)7r°(p2)|7*(9i, Ai)7*(<72, A 2 )), (5.11) 

hence 


i = 1 , 
i yf 1 , 


i = 3, 
i yf 3. 


(5.9) 


(5.10) 




1 cr^js) 

327r2 2 


dn'i 


{pi,P2,qi) M^+-*(pi,732, -<73) 

+ (731,732, gi)VhJ;r*( 73 i, 732 , -93)^ , 


(5.12) 


where the subscripts {H—, 00} denote the pion charges. The analytic continuation of the unitarity relation can 
be obtained if the 7 * 7 * —t tttt matrix element is expressed in terms of the fixed-s dispersion relation (2.42) 
for its scalar functions: 
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Ui — u 


tl — t 


dui 


’AMI 


Ui — U 


(5.13) 


Note that does not contain any pole terms because the photon does not couple to two neutral pions due 
to angular momentum conservation and Bose symmetry. 

If we pick the contribution of the pole terms on both sides of the cut, we single out box topologies: 




\<7 


_ 1 <T71-(s) 

box 327r2 2 


dVt': 
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iJ = l,A 



t' -Ml 


(5.14) 


where the primed variables belong to the sub-process on the left-hand side and the double-primed variables to 
the sub-process on the right-hand side of the cut. 

We could now apply a tensor reduction to obtain phase-space integrals for the reduced scalar quantities. The 
projection on the scalar functions 11 ^ would allow us to identify the discontinuities due to box structures. 
The reduced scalar integrals could then be transformed into another dispersive integral. Together with the 
dispersion integral ds' of the primary cut, this produces the double-spectral representation. The case of the 
simplest scalar phase-space integral is explained in [33]. 

Unfortunately, the tensor reduction is enormously complex because it involves the inversion of a 138 x 138 
square matrix. We will therefore follow a different strategy for the box topologies. We note that the non-zero 
pole residues pl'.^~ contain two electromagnetic pion form factors for the off-shell photons. These form factors 
can be factored out and multiply then the discontinuity that would be obtained by applying Cutkosky’s rules [63] 
to the sQED pion loop calculation. This becomes clear from the relation between the pole terms and the sQED 
Born terms as discussed Sect. 2.6. Therefore, the box contribution is nothing else but the sQED contribution 
multiplied by a vector form factor (qf) for each of the off-shell photons. As in the case of the sub-process, 
the difference between unitarity diagrams and Feynman diagrams is absolutely crucial: the sQED contribution 
consists of boxes, triangles, and bulb Feynman diagrams, but corresponds to the pure box topology in terms of 
unitarity. In Sect. 5.3, we will prove that this identification is correct and unique. 

Finally, there are the contributions with discontinuities either in one or in both of the sub-processes. They 
contain for example resonance contributions in the sub-process, but also two-pion rescattering effects in the 
direct channel. 
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Figure 6: Unitarity diagram representing the pion-pole contribution in one channel. 


5.2 Symmetrization and classification into topologies 

In the previous subsection, we have explained how the double-spectral representation can be derived from a 
fixed-t dispersion relation by taking the analytic continuation in t, which is defined by the unitarity relation. 
In the s-channel unitarity relation, a fixed-s dispersion relation of the sub-process is inserted (in the unitarity 
relation for the u-channel contribution, the variable u is kept fixed, which, however, plays again the role of s in 
the sub-process). Of course, one could have started with a fixed-u or fixed-s dispersion relation in the first place. 
The requirement that this lead to the same result allows us to identify a symmetric representation, which treats 
the Mandelstam variables on an equal footing and therefore implements crossing symmetry. In this symmetric 
representation, we classify the different contributions in terms of topologies. Note that in the case of HLbL, we 
get the two other possibilities (i.e. taking fixed-M and fixed-s dispersion relations as the starting point) for free, 
because we consider a totally crossing symmetric process. 

If we compare the three different representations of unsubtracted dispersion relations (5.4), we immediately 
see that the contributions in one representation are either contained explicitly in the other representations, 
can be understood as part of the respective constant c\, or correspond to a contribution of neglected higher 
intermediate states, as we are now going to discuss. 

5.2.1 Pion-pole contribution 

First, we consider the pion-pole topology (see Fig. 6). The fixed-t dispersion relation contains the poles in the s- 
and u-channel explicitly. Analogously, the fixed-s (fixed-rt) dispersion relation contains the poles in the t- and 
M-channel (s- and t-channel). The t-channel pole contribution, which is not explicit in the fixed-t representation, 
can be identified with c* as it vanishes in the limit t —> oo. Hence, the total pion-pole contribution is just given 

by 


n 


TT^-pole 


(s,t,M) = 


Pi\s 


Pi\t 


Pi;t> 


s-M^ t-M^ u-Ml' 
where the pole residues are products of pion transition form factors: 


Pi,s — ^il “^7r07*7* (^l 7 92)'^7r°7*7* (93 ) 94 ) > 
Pi,t ^i2 (Qi 7 Qs ) (^2 ? ^4)7 

Pi,u — ^ 7r07*7*(Qi7<?4).F' (^Q 2 y Qs') y 


(5.15) 


(5.16) 


where Sij is the Kronecker delta. 

Since the Lorentz structures of the usual pion-pole contribution, see e.g. [19], coincide with the first three 
BTT structures, this proves that in a dispersive approach the pion pole is unambiguously defined as given in [26] 
(and already in [19]), in particular with the transition form factor defined for an on-shell pion. 

5.2.2 Box contribution 

For the contribution of box topologies (see Fig. 7), Mandelstam diagrams prove very useful for the discussion 
of double-spectral regions. In Fig. 8, such a diagram is shown for the case = 0.5M^. A dashed line indicates 
a line of fixed t, used for writing the fixed-t dispersion relation. The two cuts are highlighted in gray. The 
discontinuities along these cuts can be written again as a dispersive integral over double-spectral functions. The 
three regions of non-vanishing double-spectral functions are labeled in Fig. 8 by pst, Psuy and ptu- 
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(a) (b) (c) 

Figure 7: Unitarity diagrams representing the box contributions. 


The s-channel cut receives contributions from the double-spectral regions pgt and according to the 
unitarity diagrams 7a and 7b, where first the vertical cut, then the horizontal cut is applied.The n-channel 
cut receives contributions again from psu and from ptu, according to the unitarity diagrams 7b and 7c. In 
diagram 7b, the horizontal cut is now applied first. 

Hence, the fixed-t dispersion relation leads to a priori four double-spectral integrals: one for each of the 
regions pst and pm and two for the region psu- However, it turns out that the sum of the two double-spectral 
integrals for the region psu equals the crossed version of one of the other double-spectral integrals. This is illus¬ 
trated for the example of a simple scalar box diagram in App. G.3. Therefore, the box contributions constructed 
from a fixed-t dispersion relation are already crossing symmetric and identical to the box contributions that are 
obtained from a fixed-s or fixed-u dispersion relation. 

This discussion already anticipates the idea of the uniqueness proof in Sect. 5.3: isolating the pion-pole 
contribution to the LHC in tttt intermediate states, the double-spectral functions become identical to the ones 
of the box diagrams in the sQED loop calculation. What needs to be shown is that the Lorentz structures in the 
BTT set not only produce the right Lorentz structures for the sQED box diagrams, but also that the triangle- 
and bulb-type sQED diagrams are correctly reproduced. Moreover, this identification has to be proven for the 
full set of BTT functions, including those which cannot immediately be identified with basis functions due to 
the occurrence of the Tarrach amplitudes. 

5.2.3 Higher intermediate states 

Topologies where we replace either one or both of the pion poles in the sub-process by a multi-particle cut are the 
last contribution with a two-pion intermediate state in the direct channel. Such contributions could be captured 
in several ways, either relying on a partial-wave picture or approximating the multi-pion states by an effective 
resonance description. While the case of S'-waves was already discussed in [26], the BTT formulation facilitates 
the generalization to D-waves, since the technical complications related to angular averages and kinematic 
singularities, the latter requiring the introduction of off-diagonal kernels, are avoided. The application to 
rescattering effects will be presented in a forthcoming publication. 

5.3 Uniqueness of the box topologies 

In this subsection, we will prove that the box topologies are equal to the FsQED contribution (sQED multiplied 
by pion vector form factors) and that this contribution is unique. The proof is based on the Mandelstam 
representation and consists of two arguments: 

• the double-spectral densities of box topologies and FsQED agree, 

• both representations fulfill an unsubtracted Mandelstam representation. 

The argument for the first point is the following: in the case of the box topologies, the discontinuities are defined 
by the unitarity relation as phase-space integrals of products of two pure pole contributions for 7 * 7 * —> tttt. In 
the case of FsQED, the discontinuities are most easily found by applying Cutkosky’s rules [63], i.e. by cutting 
the loop diagrams. We see immediately that the discontinuities in FsQED are phase-space integrals of products 
of two Born terms for 7 * 7 * —> tttt. In Sect. 2.6, we have already shown that the Born contribution is the same as 
the pure pole, hence the discontinuities (and therefore also the double-spectral densities) of the box topologies 
and FsQED are equal. 

^■^In fact, each of the shown diagrams corresponds to two topologies because the pion is charged and its line has a direction. 
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Figure 8: Mandelstam diagram for HLbL scattering for the case qf = 0.5M^ with double-spectral regions for box 
topologies. The dashed line marks a line of fixed t with its s- and rt-channel cuts highlighted in gray. 


To complete the proof it remains to be shown that the scalar functions in FsQED fulfill a Mandelstam 
representation (the scalar functions of the box topologies are defined by the Mandelstam representation). Apart 
from the first six scalar functions Hi,..., Ilg, which are unaffected by the Tarrach amplitudes, the scalar BTT 
functions Hi are defined only up to the ambiguity (3.22), so that we cannot compare the representations for 
these functions immediately. Therefore, we will derive the double-spectral representation of the basis functions 
Hi that is implied by an unsubtracted Mandelstam representation for the scalar functions Hi and show that the 
FsQED basis functions fulfill a double-spectral representation of this form. This will then complete the proof 
of the uniqueness and the equality of the box topologies and the FsQED contribution. 

First, we calculate the fully off-shell sQED loop contribution. It consists of six box diagrams, twelve triangles, 
and three bulb diagrams: 


ie'^U 


sQED 




The three bulb diagrams are given by 



M^) + g>^^g‘'^Boit, Ml) 


+ g^'^g-'^Bo{u,MlMl 


))■ 


(5.17) 


(5.18) 
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where the scalar loop function Bq is defined in (G.l). One of the twelve triangle diagrams can be represented 
as 



= 4^eV" [i^qU 2 Co{qi, s, qi M^, M^) - {q’i - ql^) , s, qi, M^, M^) 

+ 92 (9f - 92^) C2{ql s, ql M^) 

- 2g''"Coo(9?, s, ql, Ml Ml Ml) + 2 ^92^ + qUl) Ci2{qls, ql Ml Ml Ml 


- 2q'^q’(Cilqt, s, ql Ml Ml Ml - 2 q^q-C 22 iqt, s, qI, Ml Ml Ml 


(5.19) 


where the scalar loop function Cq is defined in (G.4). For the tensor coefficient functions Cj, Cij, we use the 
convention of [64] up to a normalization factor: Idn^Co = etc. The remaining eleven triangles follow by 
crossing the external momenta. Similarly, the fully off-shell formula for the box diagrams can also be expressed 
in terms of the scalar four-point function Dq and the tensor coefficient functions Di, Dij, Dijk, Dijki, but the 
result is rather long and need not be reproduced here. 

In this way, we obtain a representation in terms of scalar loop functions Bq, Cq, Dq, and the pertinent tensor 
coefficient functions. Next, we project this expression onto the HLbL tensor basis (3.18), see App. G, in order to 
identify the sQED contribution to the basis functions fti. Using FeynGalc [65], we perform a Passarino-Veltman 
reduction [66, 67] of the tensor coefficient functions, so that the basis functions are given by linear combinations 
of the scalar loop functions: 


nf = p^ 


+ CLi^oiMl 

blBoiqlMl Ml + b^Boiql Ml Ml + bfBo{ql Ml Ml + bfBoiql Ml Ml 
btBol, Ml Ml + blBoit, Ml Ml + b^B^u, Ml Ml 

clClql ql s. Ml Ml Ml + c^C.iql ql t. Ml Ml Ml + c^C.iql ql u, Ml Ml Ml 
clCoiql ql s, Ml Ml Ml + cfClql ql t, Ml Ml Ml + cfCoiql ql u, Ml Ml Ml 
dfDoiql qlql ql s, t, Ml Ml Ml Ml 
dlDoiqlql ql ql s, u, Ml Ml Ml Ml 

dlDoiql qlql ql t, u, Ml Ml Ml Ml, (5.20) 


where the coefficients pi, ai, b), c), d) are meromorphic functions of the virtualities qf and the Mandelstam 
variables s, t, u. The Passarino-Veltman reduction results unfortunately in extremely long expressions for 
the coefficients, which, in addition, contain many kinematic singularities. However, at the kinematic points 
of these singularities the scalar loop functions fulfill linear relations which cancel most of the singularities. 
By evaluating the expressions for the functions numerically, we have checked that the only surviving 

kinematic singularities are exactly the ones required by the projection of the 54 scalar functions H^ onto the 
basis (B.2). 

As already mentioned, the verification of the Mandelstam representation is simplest for the functions 
Hi,... ,n6 , because they are identical to the basis functions Fti,... ,Il6, i.e. their projection does not involve 
kinematic singularities. We have to show that these functions fulfill a Mandelstam representation of the form 
(suppressing the ie) 


n; 


sQED 




{s,t,u) = — 


ds' 


/4M2 
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r ., 

^ {s'-s){u'-u) 


4 dt' 


du' 




7^^ 94M2 Ju+{t') {f -t){u' -uY 

where the borders of the double-spectral regions ,u^ are defined in (G.12). 


(5.21) 


31 



As the scalar four-point functions Dq are the only ones in (5.20) with a double-spectral region, we can 
immediately identify the double-spectral densities: 

= df{s,t)pois,t;ql,qj,ql,ql), ') 

Pl^u^is,u) = dt^{s,u)pois,u]ql,q^,ql,ql), > t 6}, 

pt^u^it,u)=dTit,u)poit,U-,qlqlql,ql) J (5.22) 

where po is the double-spectral density of the Dq scalar loop function, defined in (G.15). By evaluating numer¬ 
ically the double-spectral integrals (5.21) at some random kinematic points below the appearance of anomalous 
thresholds, we have checked that the Mandelstam representation for the functions ..., indeed 

agrees with the expression in terms of scalar loop functions (5.20). Even though showing the identity of the two 
expressions algebraically would have been desirable, already the sheer size of these expressions has prevented 
an analytic comparison. However, the numerical check does not leave any room open for differences. 

This result implies that the double-spectral densities of the sQED box diagrams (5.22) have the correct form 
that, when inserted into the Mandelstam representation, the non-residue pieces can be separated in such a way 
that all triangle and bulb loop functions in (5.20) are reproduced with the correct coefficients. The explicit 
verification of this property is the crucial point in the argument. Since the double-spectral densities of box 
topologies and FsQED agree, this proves the claim for i = 1,..., 6 . 

The remaining basis functions exhibit kinematic singularities. In order to work out their double-dispersive 
representation, we can limit ourselves to the seven remaining representatives of (3.19). All other basis functions 
are related by crossing symmetry (B.3). 

We write the basis functions of interest as 


- 4^3 • q4ll'i = 2(s - Sh)n 7 - {u- Ua){t - t&)n 3 i, 

- 2^3 • q4flg = (s- Sb)n 9 + (u- Ub)n 22 , 

- 4^3 • 941119 = 2 (s - Sh)ni 9 -k (U - Ua)(u - Ub)n3i, 

-4qi ■ 92® • <741421 = (s - Sa)(s - St,)n 21 - (u - Ua)(u - Ub)n22, 

2qi ■ 921436 = (s — Sa)n43 — (u — '«h)n37, 

— 293 • 941439 = (s — St,)n 49 — (s — Sa)n 54 , 

—293 • 941440 = (s — St,)n 5 o — (i — tb)n54, (5.23) 

where Sa '■= qf + 92 , := <73 + <? 4 , Ua '■= 92 + < 73 , Ub ■= < 7 i + <? 4 , h := 9 ^ + 94 . Next, we introduce subtractions 

for the scalar functions Hp each function is subtracted once or twice at the subtraction point that appears in 
the coefficient of the other function H^ in the above linear combinations. For instance, in the equation for n 9 , 
we subtract 09 at u = Ub and II 22 at s = Sb- This leads to representations of the basis functions Hi where both 

double-spectral contributions are multiplied by a common prefactor (e.g. (s — Sb){u — Ub) in the case of n 9 ), 

while only the subtraction terms, which are single-dispersion integrals, have different prefactors. The results 
for these double-spectral representations of the basis functions fti are given in (F.2) and (F.5) in App. F.l. 

The important point is that in these double-spectral representations both the discontinuities of the single¬ 
dispersion integrals and the combined double-spectral densities pi are unambiguously defined. We do not 
attempt to split the double-spectral densities into two contributions from the H^: this splitting involves exactly 
the ambiguity (3.22). For the sQED contribution, the discontinuities and double-spectral densities can be 
extracted from the loop representation of the basis functions (5.20) as follows: 

• The discontinuities are extracted by taking the appropriate limit of the basis functions H^. For instance, 
the discontinuity of the subtraction term of Hg is obtained by considering the limit lim (—293 • 94 n 9 ). 

U—¥Uh 

• The combined double-spectral densities pi are (up to a kinematic factor) just the double-spectral densities 
of the basis functions 14^. Therefore, they can be obtained again from the coefficients of the Dg functions 
in the loop representation, in analogy to the case of Hi,..., He. 

After having identified all the discontinuities and double-spectral densities in (F.2) and (F.5), we have checked 
numerically for random kinematic points (below the appearance of anomalous thresholds) that the dispersive 
representations of the functions 14^ agrees with the loop representation. It turns out that = 0, 

hence we can set 

= ... = = 0, (5.24) 
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which also fixes the redundancy (3.23). 

This completes our proof of the uniqueness of the pion-box contribution. Cutkosky’s rules tell us that the 
discontinuities of the FsQED contribution are the same as the ones of the pion-box topologies in the sense 
of unitarity. The FsQED contribution fulfills the same double-spectral representation as the pure pion-box 
topologies. Therefore, the two representations are the same. Unitarity and Mandelstam analyticity define the 
pion-box contribution, i.e. tttt intermediate states with a pion-pole LHC, in a unique way. 

Finally, we stress that these calculations also provide a strong test of the Lorentz decomposition (3.13). 
Apart from the function n 49 , which does not receive a contribution from the pion loop, all scalar functions have 
been shown to be free of kinematics, so that the sQED pion-loop amplitude behaves exactly as expected from 
the general BTT formalism. 


5.4 Contribution to {g — 2)^ 

In this subsection, we insert our dispersive representation of the scalar functions into the master formula (4.28) 
to obtain the contribution to a^. 


5.4.1 Pion-pole contribution 


With (5.15) and using the master formula (4.28), we recover the well-known result for the pion-pole contribution 
to [19]: 

rt 3 nOO pOO pi 

I dQi dQ2 J ^drVl-r^QlQl 

X (Ti(gi,g2,r)nf-p°''’(gi,Q2,r) + T2(Qi,Q2,T)nf-p°“^(gi,Q2,T)), ( 5 . 25 ) 

with 
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r®-pole _ ( QlJ Q2)‘^7r®7*7* ( Q 37 O) 


gi + M2 
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r°-pole _ .77^07*7* ( gii { g 27 

' “ Ql + M 2 ’ 

where g| = gf + 2 Q 1 Q 2 T + g| and the integral kernels Ti are given in App. E.2. 


(5.26) 


5.4.2 Pion-box contribution 

The single-integral discontinuities and the double-spectral densities in the dispersive representations of the basis 
functions (F.2) and (F.5) are quantities that can be extracted directly from the projected basis functions 11^. 
To the contrary, the separation of the double-spectral densities pi into the two contributions from the different 
scalar functions H^ is not unambiguously possible, which reflects just the redundancy (3.22). However, such a 
separation is not necessary: for the calculation of a^, we need the scalar functions H^ only in the limit fc —> 0. In 
this limit, all the scalar functions H^ appearing in the master formula (D.2) can be expressed in terms of single¬ 
dispersion integrals, where the discontinuities are directly related to the basis functions 11^. All the subtracted 
double-spectral integrals, which are not unambiguously defined, drop out in the limit A: —)■ 0. 

We stress that the ambiguities in the definition of the BTT functions H^ are related to the fact that these 
functions are not observables. In the calculation of any physical quantity, such as helicity amplitudes or the 
HLbL contribution to a^, any ambiguity due to the redundancy in the BTT set has to drop out. 

The pion-box contribution to is therefore given by 

ey 3 poo poo pi 

^.-box ^ ^FsQED = —j^ dQi dQ2 j ^ drVl - T^QlQl i-Ql)F^{-Ql)F^i-Ql) 

X VT,(gi,g2,T)nf^°(gi,g2,T), 

(5.27) 

where the functions H^ are defined in (D.2). They are linear combinations of the scalar functions H^ in the 
limit k —>■ 0. The required functions H^ can be obtained in this limit from the basis functions H^, see (F.8) in 
App. F.2. Their explicit dispersive representation is given in (F.9). In fact, this representation is applicable not 
only to the pion-box contribution, but whenever a dispersive representation for the H^ can be constructed. In 
the special case of the pion box, one has in addition = 0, hence = 0. 
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6 Conclusion and outlook 


We have derived a decomposition of the HLbL tensor into scalar functions (3.13), following the general recipe by 
Bardeen, Tung, and Tarrach, which allowed us to derive a master formula (4.28) for the HLbL contribution to 
the anomalous magnetic moment of the muon {g — 2)^ in terms of scalar functions that by construction are free 
of kinematic singularities and zeros. In particular, these scalar functions fulfill a Mandelstam double-spectral 
representation. In the present work, we have considered only the lowest-lying intermediate states, pion-pole and 
pion-box topologies. These two contributions can be studied without any approximation. The central result is 
that the pion pole corresponds exactly to the contribution of one-pion intermediate states if interpreted with 
an on-shell pion transition form factor, and the pion-box, i.e. tttt intermediate states with a pion-pole LHC, to 
the sQED pion loop dressed with pion vector form factors. These results confirm both the definition of the pion 
pole and prove explicitly the validity of the separation of Born and non-Born tttt contributions given in [26]. 

Moreover, while the formulation in terms of helicity amplitudes adopted in [26] required the introduction of 
off-diagonal kernel functions to account for kinematic singularities, the BTT version avoids these complications 
and thus facilitates the generalization to partial waves beyond S'-waves. In order to exploit unitarity in a 
partial-wave picture, the partial-wave unitarity relations analyzed in the helicity basis need to be transformed 
to the BTT basis, the formalism for which we provided in the present paper. 

This treatment is based on fundamental principles of particle physics: gauge invariance and crossing sym¬ 
metry are already implemented in the decomposition of the HLbL tensor into scalar functions. The dispersive 
description uses analyticity and unitarity to establish a relation between the HLbL contribution and different 
on-shell quantities. These on-shell quantities are in principle either experimentally accessible or can be recon¬ 
structed from data with dispersive methods. For the two cases covered in this article, the corresponding input 
is parametrized by the pion transition form factor in case of the pion pole, and by the pion vector form 

factor for the pion box, both form factors being needed for negative virtualities of the off-shell photons. For 
general tttt intermediate states also information on the partial waves for 7 * 7 * —)■ tttt is required. An overview 
over which processes can help in the dispersive reconstruction of the pion transition form factor [68-71] or 
7 * 7 * —> TTTT [26, 36, 40, 56, 72] is provided in [27]. 

An extension of the presented dispersive treatment is possible within certain limits. It is straightforward 
to include higher pseudoscalar poles, i.e. the ry and rj' mesons by just adding their contribution in complete 
analogy to the 7 r° pole. The input quantities will be the transition form factors of these heavier pseudoscalars, 
see [73-75] for work towards their dispersive calculation. Although these mesons are unstable in QCD, their 
decay width is certainly small enough to justify the treatment as a pure pole. A bit more difficult is the inclusion 
of higher two-particle intermediate states: an extension to KK intermediate states is still straightforward, but 
the experimental input will be less accurate. In general, the approach applies to the low and intermediate 
energies < 1.5 GeV, where a few channels dominate. All model calculations done so far indeed support the 
assumption that the lowest-lying singularities govern the HLbL tensor. 

A careful numerical analysis of the formalism presented in this paper is in progress. It will reveal the relative 
importance of the different contributions and should allow us to identify which input quantities will have the 
largest impact concerning the reduction of the hadronic uncertainty in In this way, we are confident that 

the presented treatment of HLbL scattering shows a path towards a data-driven and thus less model-dependent 
evaluation of (g — 2 )^. 
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A Partial-wave representation for 7*7* —)■ tttt 

Here, we give the expressions for the scalar basis functions of the sub-process 7 * 7 * —> tttt, defined in (2.26), in 
terms of the helicity partial waves including I?-waves: 


Ai — 


Ai2(s) 


1291^2 (^*o.5(s) + - 1 )^ 2 , 5 ( 5 )^ -(s-q^- qj) ^ho,i(s) -f ^(30^ - l)h 2 ,i(s)^ 

5v^ {{{ql - qlf - s{ql + ql)) - s{s - qf - ql)) 


4s 


-^ 2 , 2 ( 3 ) 


^ (s(9l + 92 ) - (9l - ql f) Z^h2,3{s) + ^(9l - 92 )(s - 9? - gi)2^/l2.4(s) 


Ao — 


Ai2(s) 


2 yho,iis) + -(3z^ - l)/i 2 .i(s)j - {s-ql- ql) (^/io,5(s) + -(32:^ - 1 )^ 2 , 5 ( 5 ) 

, ,, , 2 ^^ t \ , K AT 2 u ( \ ( 2 ( 9 ? - <?i)^ - Ai 2 (s)) ^ 2 , 4 ( 5 ) 

+ 5 \/ 7 1 + Z^)h2,2{s) + 5V3sZ^h2,3[s) -^^- 2 ^- 2 

V 2 2 vs of - 05 


91 - 9 i I ’ 


^3 = - 

A4 = — 

= 


5v^ 


(s-4M2)Ai2(s) 


^ 2 , 2 ( 5 ) + \/ 7 / 12 . 3 ( 5 ) - \/ 7(91 + 92 ) „2^’"^^„2 ; ’ 


91 - 91 

loVe 


. 2 2 

91 - 92 


(s-4M2)Ai2(5) 


(^/l 2 . 2 (s) + y|(s- 9 ?- 92 ) 


/l2.4(s) 
2 2 
91 - 92 


(A.l) 


We note that under q\ o ql, /i 2,4 changes sign, while all the other S- and I?-waves are invariant. Therefore, 
the apparent singularities in ql — ql in (A.l) are canceled by the corresponding kinematic zero in /i 2 . 4 . 


B Tensor basis and crossing relations 

B.l Crossing relations between Lorentz structures 

The redundant set of generating Lorentz structures of the HLbL tensor consists of 54 elements. While seven 
structures are defined in (3.14), all the remaining structures are given by crossed versions thereof: 

rp^uXa 
rp^uXa 
rp^uXa 
mUi^Xa 


pfiuXcr 

^20 

p ^ vX <7 

^23 

pflvXcT 

^26 

p ( lvX <7 

^29 

p^vX<7 

^32 


uXg 


‘-35 


Ci4[rr^1, 

rpfluXcr 

-^3 

=ci3[Tr^"], 



ci4[Tr^i, 

rpfluXcr 

-^6 

=ci3[Tr^"]: 



Ci2[Tr^1, 

rpfluXcT 

^9 

= Ci3[C23[T7^"^"]], 

rpIluXa 

-^10 

= C23[Tr^1, 

c24[Tr^i, 

rpfluXcr 

-^12 

= Ci4[C24[T7^"^"]], 

rp(lvX(7 

^13 

= Ci3[Tr^"], 

C23[Cl3[r7''"^"]], 

rpfluXcr 

^15 

= Ci4[Tr^"], 

rplAUXcT 

^16 

= C24[Ci4[r7^"^": 

C24[Cl3[r7''"^"]], 

rpfluXcr 

^IS 

= C23[Ci4[Tr^"]], 



C34[r|(r^"], 

rpfluXcr 

-^21 

= C23[Ti7^n, 

rpIluXa 

-^22 

P \P \rjilll'X<7 

— <-24 [<-23 [-^19 

C23[C24[rr9"^"]], 

rpfluXcr 

= C24[Ti7^4, 

rpIluXa 

^23 

p \pt \rjilll'X<7 
— <-23 [<-13 [4 19 


rpfluXcr 

^27 

= Ci4[Ti7^"], 

rpIluXa 

-^28 

= C24[Ci4['ri7^" 

C24[Ci3[Ti7^"]], 

rpfluXcr 

^30 

= C34[C24[Ci3[Ti7^1]], 



C24[Ci3[T37^"]], 

rpfluXcr 

^33 

= C23[T3T^n, 

rpIluXa 

^34 

= Ci3[T37^"], 

C24[T37^4, 

rpfluXcr 

^36 

= Ci4[T3T^n, 



<-34 [<-14 [-^37 JJ) 

rpfluXcT 

^39 

p rrplAl'' X(T '\ 

— <-14[-^37 \i 

rpIluXa 

-^49 

P \P \rjilll'X<7' 
— <-12 [<-14 [4 37 

C23[Cl2[r3T^"]], 

rpfluXcT 

-^42 

= Ci2[C24[T37^"]], 

rpIluXa 

^43 

= C24[T37^4, 

Cl2[C23[r37^"]], 

rpfluXcr 

-^45 

= C23[T37^"], 

rpIluXa 

-^46 

= Ci4[C23['r37^": 
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T \[IV\<7 /-t \rf\fiu\<7'\ 

52 -^ 13 [-t 49 J 


rplll>X<7 

= C24[Ci3[T37^"]], 

rpfluXcr 

-^48 — 

rp(ll>\(7 

-^50 

= Ci2[C24[IT9"^"]]. 

rpfluXcr 
^51 — 


= Ci2[Ci3[T';,"^"]], 

rpfluXcr 
-^54 — 


Ci2[T37^n, 

C2:i[CiA[T^t"]]- 


(B.l) 


B.2 Basis coefficient functions 

The relation between the 43 basis coefficient functions 11^ and the 54 scalar functions of the redundant set can 


be obtained by projecting (3.13) on the basis. 

The result defines the exact form of the kinematic singularities 

in qi ■ q 2 and 53 • (74 in the functions 11 ^: 




Hi = Hi, 

n 2 

= n 2 , 

n 3 = 113 , 

n 4 = 114 , 

fis 

= IIs, 

fie = Ilg, 

Hr = Hr - 

fis 

jr 9i-939i-94„ 

= iis iisi. 


93 • 94 


93 • 94 


fig = Hg + ■ '^"n 22 , 

fllO 

_ rr 1 92 ■ 93 TT 

— AAlO + AI 22 , 


93 • 94 


qi • q 2 


fin = Hn - ^^n24, 

ni 2 

= ni 2 - ^^n24, 


9i • 92 


93 • 94 


ni3 = ni3 + ^^n26, 

ni4 

— 1114 + II 26 , 


93 • 94 


9i • 92 


fr _ TT 92 ■ 93 TT 

file 

— TT 9i ■ 94 TT 

— 1 I 16 1 I 28 ; 


93 • 94 


9i • 92 


filT = ni7 - • g4g2 • 94 

fiis 

„ 9i • 9392 • 93„ 

— II 18 1132) 


9i • 92 


9i • 92 


f, TT , 9i-9492 - 93„ 

ill 9 — liig + 1131 , 

n2o 

„ 91-9392-94 

— 1120 + 1131 ) 


93 • 94 


93 • 94 


n21 = n2i - • g4g2 • 93 

n22 

jr 9i-9392 - 94„ 

— II 23 1124) 


9i • 9293 • 94 


9i • 9293 • 94 


fr „ 91-9392 - 94 

^^23 — J^i25 J^J^26 7 

n24 

9i • 9492 • 93 JJ 

— 1127 1128 ) 


9i • 9293 • 94 


9i • 9293 • 94 


fr „ 9i-9492 - 93„ 

1125 — AI 29 AA 327 

n26 

„ 9i-9392 - 94„ 

— AA 30 A 132 , 


qi • q2 


qi • 92 


fi27 = n33 + n22, 

n28 

= n34 + n26, 


9i • 9293 • 94 


9i • 9293 • 94 


n29 = Has-n24, 

fiso 

n 91 • 93 „ 

— AA 36 


9i • 9293 • 94 


9i • 9293 • 94 


n3i = Has + ^^n47, 

n32 

= n39 - ^^n46, 

n33 = n4o - ^^n46) 

93 • 94 


93 • 94 

93 • 94 

9i • 93 

n34 = n4i + ^^^n47, 

flss 

= n42 + ^^n48, 

636 —1143 + 1137 , 

93 • 94 


9i • 92 

9i • 92 

n37 = n44 '^"■®n48, 

fiss 

= n45 - ^^n37, 


9l • 92 


9i • 92 


£[39 — 1149 + 1154 , 

n4o 

= Hso - ^^n54, 

n4i = iisi + 1154 , 

93 • 94 


93 • 94 

93 • 94 

n42 = n52 - ^^n54, 

n43 

= n53 + ^^n54. 

(B.2) 

93 • 94 


93 • 94 

Crossing symmetry for the functions 11^ 

implies the following relations between the basis functions 11^: 

^2 = Ci4[ni]) 

63 

= C 13 P 1 ], 


fis = Ci4[n4], 

fie 

= C 13 P 4 ], 


fis = Ci2[n7]) 

ni7 

= C24[Ci3[n7]], 

fflS = C23[Cl4[n7]]) 

filO = C24[Cl3[n9]], 

fill 

= C34[C24[Ci3[n9]]], 

ni2 = C34[n9], 

ftis = Ci2[fi9], 

ni4 

= Ci2[C24[Ci3[fi9]]]) 

filS = C34[Ci2[n9]]) 

file = C23[Cl4[n9]]) 
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n27 = 


Ilao = — 


C 34 P 21 ], 

1123 

C34[ni9], 

n25 

^23^34 pig]]. 

n28 

9i • 9293 • 94 


C24[Cl4pi9]], 


9i • 9293 • 94 


C24[Cl3p36]], 

£[32 

Cl2[C24[Cl3p36]]], 

£[35 

C 34 Pse], 


Cl2p4o], 

£[42 


<?1 • 94 


9i • 9293 • 94 


-Ci3[n 


19j) 


1124 — C34[Ci2[n2i]], 
n26 = C34[C23[Ci4[nig]]], 
92 • 93 


1129 = — ■ 


9l • 9293 • 94 


C 24 P 


19j) 


£[33 = C23[Cl4[n36]], 
£[37 = C34[Ci2[n36]], 

£[43 = C34[Ci2[n4o]]- 


(B.3) 


1142 = C 34 P 40 ], 

C Projection of the scalar functions 

Given any representation of the HLbL tensor the following procedure allows the identification of the 

basis coefficients 11^ in (3.18): 

• write the tensor in terms of the 138 elementary structures (3.8) and identify the scalar coefficients 

• take the subset consisting of the following 43 scalar coefficients: 

{SJ = {Sp = 1,2,3,4,5,7,8,16,17,19,20,25,26,28,29, 

34,35,37,38,43,44,46,47,53,54,56,57, 

94,95,97,98,103,104,106,107, 

121,122,124,125,130,131,133,134}, 


(C.l) 


perform a basis change according to 


43 


Uj — ^ 


where P is an invertible 43 x 43 matrix, defined below. 

The matrix describing the change of basis is sparse and has the following non-zero entries: 


912934 


P2,3 = 


^3,2 = 


714,4 = — 


76,4 = — 


7*8,2 = 


1 

7*2,4 = 

913924 1 1 

7^2.5 = 

1 

7*2,19 = 

1 

7*2,25 = 

1 

912934 ’ 

9i29i4 912934 ’ 

912934 ’ 

9?2934 ’ 

912914 ’ 

1 

7*3,4 = 

914923 1 1 

7*3,6 = 

1 

7*3,20 = 

1 

7*3,26 = 

1 

912934 ’ 

9?29i4 912934 ’ 

912934 ’ 

9?2934 ’ 

912914 ’ 


9l39l4 

912 9 I 4 ’ 
914923 

912914 ’ 

913 


7*4,7 = 


7*6.20 — 


912934 

1 


7*5,4 = — 


913924 
912 9 I 4 ’ 


7*5,19 — 


912934 

1 


7*8,17 = — 


2912934 

914 


7*8,4 = 


912934 
9l3 914923913 

2912934 9?29i4 


„ 923924 „ 

-^7.4 =-^7,8 = 

9129 I 4 912934 


7*8,6 = 


9l3 


2912934 


7*8,7 = — 


923 

9 ^ 2934 ’ 


912 9 I 4 ’ 


7*8,32 — — 


1 


7*8.35 = —: 


1 


7*9,3 = 


923 


7*9.19 = — 


2912934 

923 


7*9,4 = 


2912934 ’ ’ 2912534 ’ 

923 , 913924923 „ 923 


7*8,40 — 


2912934 ’ 


7*9,5 = 


2912934 


7*9,17 = — 


924 


912914 ’ 


7*10.3 — 


9^2934 ' 
914 

2912934 ’ 


2912934 9?29i4 

7*9,33 = X-, 7*9^36 = —--, 7*9^41 = -- 

2912934 2(712 9 3 4 2(712(734 

D 9l4 , 913924914 n 9l4 n 924 

.C^10,4 — 


2912934 9?29i4 


7*10.5 = ^^^, 7*19.7 = - 


X-’ ^10.7 

2912934 


912934 


(C.2) 
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-Pl0,25 

-Pll,3 

Pll,2i 
Pl2,2 
Pl2,26 

Pl3,2 = 


Pl3,23 — -- 


gl4 

9l2<?i4 ’ 

<724 

<712(734 ’ 

1 

<7l2<734 ’ 
<7l3 

2<7l2<734 ’ 
gl3 

9l2<7|4 ’ 

g23 

<7l2<734 ’ 

1 


Ao,31 - 
Pll,4 = 
Ai,25 = 
^12,4 = 
Pl2,32 = 
Pl3,4 = 


1 


-Pin.37 — — 


1 


o ) -* 10,37 

2<712<734 

<724 


gl3<724 
912934 
924 

912914’ 

9l3 

2(712934 

1 


2^12934 ’ 
924 


Pl0,42 — —: 


1 


912934 


Pll,5 — 


912934 


Pll,19 = — 


2<712934 
924 


9?2934’ 


914923913 
2 2 ’ 
912934 

Pi 2.3,1 = 


p _ 913 P _ 

-^12,6 — i;-' 4^12.7 — — 


O -’ ^12.7 

2912934 


923 

9?2934 ’ 


1 


r, ) 12,35 

2gi2^34 

Q23 


914923 

9?29|4 


2912934 ’ 

923 


Pl2,40 — 


1 


2912934 


912934 


Pl3,6 = 


912934 


Pl3,20 — ~ 


923 

9?2934 ’ 


Pl4,3 

Pl4,18 

Pl5,2 


912934 

914 

2912934 ’ 

9l3 


912914 ’ 


924 


Pl5,20 — — 


2912934 

924 


9?2934 ' 


Plfi.2 = 


9l3 


16,2 

Pl6,20 

Pl7,3 
Pl7,17 = 

Pl8,2 = 
Pl8,22 = 

Pl9,2 = 
Pl9,26 = 
P20,3 = 


2912934 

9l3 

9^2934 ’ 
923 

2912934 ’ 

924 

912 9 I 4 ’ 


P20,21 — —■ 


914 

912934 ’ 

1 

912934 ’ 
924 

2912934 ’ 

924 

912914 ’ 
9l3 

912934 ’ 

1 


P21,3 — 


P 2 I .25 — 


912934 

923 

2912934 ’ 

923 

912914 ’ 


Pl3,26 — 
Pl4,4 = 
Pl4,31 
Pl5,4 = 
Pl5,34 
Pl6,4 = 
Pl6,32 
Pl7,4 = 
Pl7,33 
Pl8,4 = 

Pl8,26 = 

Pl9,4 = 
Pl9,34 = 
P20,4 = 
P20,25 
P21,4 = 


923 

912914’ 

914 


913924914 

9?29i4 


Pi 4.37 = 


2912934 

1 

= ;;-> -f'14,37 

2912934 

924 , 914923924 

2912934 


Pl4,5 = 


914 


2912934 ’ 


p _ 924 

^14,7 —- 2 -’ 

912934 


1 


2912934 


, Pl4,42 = — 


1 


2912934 


9?29i4 


p _ 924 

2912934 


Pl5,18 = — 


1 


Pl5,38 — 


1 


2912934 

913 , 914923913 


9?29i4 


Pi 6.3,5 — — : 


2912934 

1 

- 14 -, -f'16,35 

2912934 

923 , 913924923 

2912934 


2912934 

Pl6,6 = 


Pl5,43 — — 1 


1 


9l3 


2912934 ’ 


2912934 

Pl6,17 = — 


923 

912 9i4’ 


914 

9129 I 4 ’ 


1 


9?29i4 


2912934 

Pl7,5 = 


Pl6,40 — 


1 


923 


2912934 ’ 


2912934 

p _ 913 

912934 


1 


2912934 


Pl7,36 — — 


1 


923914 

9 L 9 I 4 

914 

912 9I4’ 

924 

2912934 

1 

2912934 ’ 

9249?3 
2 2 
<1 i2<iL 

9l3 


914 

912934 ’ 


Pl8,6 = 


2912934 

914 

912934 ’ 


Pl7,41 = 


1 


Pl8,20 — ~ 


2912934 

914 


9?2934 ’ 


914923924 
2 2 ’ 
912934 

Pi 9.38 = — 


p _ 924 p _ 9l4 

-f^l9,6 — , -f^l9,8 —- 2 -’ 

2912934 912934 


1 


9l3 

912934 ’ 


2912934 ’ 

9l3 


Pl9,43 — — 1 


1 


P20,5 = 


P 


2912934 

9l3 


912934 


20,19 — 


9?2934 ’ 


P21,33 — 


912914 ’ 
923 

2912934 

1 


P22,3 — 


9l4 


2912934 


P22,4 — 


2912934 

914 


913924923 
2 2 ’ 
912934 

P21,36 = 

913924914 


P21,5 — 


923 n _ 9l3 

-) -f^21,8 —- 2 -’ 

2912934 912934 


1 


2912934 


) P21,41 — 


1 


9?29i4 


P22,19 — ~ 


P23,2 — 


9l4 

9?2934 ’ 
924 

2912934 ’ 


P22,31 — ~ 


P23,4 — 


2912934 

1 

2912934 

924 914923924 

2912934 9 ? 29 i 4 


P 22 5 — 


9l4 


P22,37 — 


2912934 

^22,5 = -1 P22,18 = ~ 

2912934 

1 p ^ _^ 

2912934’ 2912934’ 


913 

912914’ 


P23,6 — 


924 p _ 9l4 

o ’ -^23,8 2 ’ 

2 ^ 12^34 ^ 12^34 
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^23,18 — 

923 

P _ 1 

P 23,38 = 

1 

-f23,43 = 

1 


912914 ’ 

-^^23,34 c* 5 

2^12^34 

2512934 ’ 

2912934 ’ 


II 

913923 

p 1 

-^25,4 = — 

913924 

-P25,25 = “ 

1 


9?2934 ’ 

-124,17 — 7 

912934 

9?2934 ’ 

912934 ’ 


-P26,4 = 

914923 

p _ 1 

-P27,4 = 

914924 

-P27,18 = 

1 


9?2934 ’ 

-( 26,26 — 7 

912934 

9^2934’ 

912934 ’ 


-P28,2 = 

9?3 

_ ®13 1 

914923913 

-P28,6 = 

9?3 

-P28,7 = — 

913923 

2912934 ’ 

4^28,4 — n- 

2gi2^34 

2 2 ’ 
952934 

2912934 ’ 

9?2934 ’ 


P2s,,n — 

^29,3 = 
-f29,17 = 
-f29,36 = 
-P30,2 = 
-^30,17 = 
-^30,35 = 

^31,3 = 
-^31,17 = 
-P32,3 = 

Pd.2,W = 

P32,37 = 

^33,3 = 
-^33,21 = 
^34,1 = 
-^34,25 = 
-^34,37 = 
^35,3 = 
-P35,16 = 
-f35,36 = 
P36,2 = 
P36,ll = 


’ll; 

913914 n 1 n 9l3 „ 9l3 

912934 912934 2(712934 2(712934 

913923 n 923924913 , 923913 „ 913923 „ 1 

2912934 9i2934 2912934 2912934 934 

913924 „ 913923 □ 923 n 9l3 

912954 912934 912934 2912934 

913 p _ 913 

2912934 ’ ’ 2912934 ’ 

913923 n 913914923 , 913923 „ 913923 „ 1 

-7 -130,4 — 2—2-1" 7;-’ -^30,6 — 7;-7 -130,13 — -7 

2912934 9i2934 2912934 2912934 934 

914923 □ 913923 □ 9l3 „ 923 

912954 912934 912934 2912934 

923 n _ 923 

o ’ -^30,40 c \ ’ 

2^^12934 2gi2g'34 

9^3 0 _ 923 , 91392491 

o ’ -^31,4 2 2 

2^12^34 2^12^34 ^12^34 

Qr>‘^QoA _ 1 _ 


<&3 , 9139249 I 3 p <l23 

-- 2—2 ’ -^31,5 — 7 ;- 

912934 2912934 


-P31,8 ■ 


912934 ’ 


2 2 

' " 1 ^yi2y34 

p 923 p 923 

-131,33 — —7;-7 -131,36- 

2912934 


2912934 

913914 n I 

7-7 -1432,7--- 

4 

9l3 


923924 „ 

-^7 -131,28 

912934 

913914 n 9149249i3 , 914913 p 913914 

7;-7 -132,4 — 2—2-1" 7;-’ ^ 32,5 — 7;- 

2912934 9i2934 2912934 2912934 

1 p 914 n 913914 n 9l3 

9l2 912934 912934 2912934 

913 p _ 9l3 

2912934’ ’ 2912934’ 

913924 n 913924 , 913924 „ 913924 , . „ 

-7 -133,4 — ^—2^ H-7 -133,5 —-H J^7 -133,19 ■ 

912934 912934 912934 912934 

924 p _ 913 p _ 913924 

-7 -1^33,24 — -7 -133,25 — -^7 

'12934 912934 

1 


924 p _ yi3 

7 -133,24 — 7 ^ 33,25 — 

912934 912934 

P _ 913914923924 p _ 1 
’ -* 34,4 — 2 2 ’ -*34,5 — ^ ’ 

952934 2 


-^34,6 — 7 7 -^34,20 


-P34,34 


913 p _ 9l^ 

2912934 2(7i2( 

924 p _ 

- 134,43 


914923 p _ 924 

- ^7 - 134,32 — 7 - 

912934 2(^(12534 

923 p _ ^ p _ 924 

2912934 2(734 2gi2934 

923924 p _ 913923924 , 923924 p _ 923924 p 

2912934 9 i 2934 2(712534 2512934 

1 P _ 923 P _ 923924 p _ 924 

9l2 912934 912934 2512534 

924 P 924 

-7 -135,41 — 7-7 

712934 2512534 

913914 


2512934 ’ 

_ 913923914 I 913914 p _ 9l39l4 

9 i 29 i 4 2512534’ 2512534’ 


913914 p _ 9139239i 
o ’ -* 36,4 2 2 

2^12^34 912^34 

Ip 913 p 913914 p 914 

-7 -1 36,22 — -7 -1 36,26 — -^7 -1 36,32 — —7- 

9l2 912934 912934 25125; 


913 

2912934 ’ 


923 

2912934 ’ 


39 



-P36,35 — 
^37,1 = 
-f37,26 = 
-^37,38 = 
^38,2 = 

Pd.&,22 = 

^39,2 = 


2^ -P37,19 


<714923 
9?2934 ’ 


914 n _ 914 

2^12934 2g'i2g34 

1 n 913914923924 „ 1 n 

;7, -f37,4 — -2—2-’ -^37,5 — 71) ^37,6 - 

2 9f2934 2 

913924 923 □ 9l4 n 924 

-^7 ^37,31 — —7;-) -f37,33 — Ti-) -n37,35 — 7^-7 

912934 2(7i 2934 2qi2934 2^12934 

913 □ 1 n 9l4 n 923 

— 7 ;- 7 -f37,39 — 7 ;- 7 -f37,41 — -> -^37,42 — ”7;- 7 

2912934 2^34 2^12934 2912934 

914923 n 914923 , 914923 „ 914923 , . „ ^1495 

-7 ^38,4 — — 2' 4-’ 4^38,6 —-H ^7 -^38,20 —- 

912934 912954 912934 912934 

923 D 9l4 n 914923 

- 7 ^38,23 —- 7 -‘ 38,26 — - ^7 

912934 912934 912954 

914924923 , 924923 „ 923924 „ 

- 2—2 - 1 " 7 ;-’ -^39,6 — 7 ;- 7 439,8 - 

912934 2^12934 2^12934 

924 p _ 923924 p _ 9 


923924 p _ 91492492 
2912934’ 9 ? 29 i 4 

1 p 

- 7 439,23 

912 

923 p 

4^9,43 — —7;- 

2912934 


-f39,14 

-^39,38 — —7;-5 

2^12^34 

Pao,z = 7 


924 p 923924 p 923 

- 7 4^39,26 — -^7 4^39,34 — 7;- 7 

912934 912934 2^12934 

923 


9?4 p 9?4 , 9139249?4 p 9m p 

-5 -^40,4 — 7;-^-2-2-’ -^40,5 — 7^-5 -^40,7 — 

?12934 2 ^ 12^34 ^ 12^34 2 ^ 12^34 

P _ 913914 p _ 1 p _ 914 p _ 914 

4^0,18 — - ^7 4 - 40,29 — - 7 4-110,31 — 7 ;- 7 4 - 110,37 — 7 ;- 

912934 912934 2912934 2912934 


-f40,42 — 


912914 

914 

2912934 ’ 


P41,3 = 

914924 

II 

913914924 

914924 

P41,5 = 

914924 

P41,13 — - 
t 

2912934 ’ 

2 2 ' 
9f2934 

2912934 ’ 

2912934 ’ 

P _ 

913924 

-P41,19 

_ 914924 

P41,24 = 

914 

-^41,31 = 

924 

-^41,18 — 

912914 ’ 

9?2934 ’ 

912934 ’ 

2912934 ’ 

p _ 

924 

P41,42 ■ 

_ 924 





41,37 — 

2912934 ’ 

2912934 ’ 





-P42,2 = 

914924 

P42,4 = 

9239249m 

924914 

P42,6 = 

914924 

P42,1‘2 = “ 

2912934 ’ 

2 2 
912934 

2912934 ’ 

2912934 ’ 

P _ 

914923 

P42,20 

_ 914924 

P42,22 = 

924 

P42,34 = 

9l4 

42,18 — 

912914 ’ 

9?2934 ’ 

912934 ’ 

2912934 ’ 

1142,38 = 

914 

P42,43 ■ 

914 





2912934 ’ 

2912934 ’ 





-P43,2 = 

124 

P43,4 = 

924 , 914923924 

-P43,6 = 

<124 

P43A = — 

2912934 ’ 

2912934 

2 2 ’ 
912934 

2912934 ’ 

-P43,18 = 

923924 

P43,30 

1 

P43,34 = 

924 



912914 ’ 

912934 ’ 

2912934 ’ 



1143,38 = 

924 

P43A3 ■ 

_ 924 





2912934 ’ 

2912934 ’ 






where Qij := qt ■ Qj. 

D Scalar functions contributing to (gf — 2)^ 

The following 19 linear combinations of 33 scalar functions Hi contribute to {g — 2 )^: 
Hi = III + 9 i • 921147, 

n2 = 112 — - (91 + 9 i ■ 92) (21147 ~ IIso — IIsi — 1154), 

As = 113 — - (91 • 92 + 9I) (21147 — IIso — IIsi + 1154), 


(C.3) 
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64 — 114 + (^l + 91 • 92) Ilig + (gi • q 2 + 92) 1120 

+ (9? + 9i • 92) (91 ■ 92 + 92) nai — 2 (21147 — IIso — IIsi) + - (91 — 9!) n54, 

ns = Ho — (71 • 921121 + 2 (91 ■ 92 + 92) (21122 — 2(7 i • 92X133 + IIso + IIsi — 1X54) — 9IXX47, 

IXe = XXg — 9i • 92XX25 + - { q \ + 9i • 92) (2XX26 — 2(7 i • 92XX34 + XX50 + XX51 + XX54) — 93XX47, 

1X7 = XX7 — XX19 — (91 • <72 + 92) 1X31, 

Hg = XXg — XX20 — (91 + 9 i • 92) XX31, 

Hg = XXg — XX22 + 9l • 92XX33, 

IXio = XXio — XX21 — (91 • 92 + 92) 1X33) 

IXii = XXii + XX47 — XX54, 

1X12 = XX13 — XX26 + 9i • 92XX34, 

1X13 = XXi 4 — XX25 — {q\ + 9i • 92) XX34, 

IXi4 = XX16 + XX47 + XX54, 

Ills = XX17 + XX47 — XX50 — XX51, 

XXie = 2 (XX39 + XX40 + XX 4 g), 

ni7 = XX42 — XX47 + - (XX50 + XX51 + XX54), 

fiig = XX43 — XX47 + - (XX50 + XX51 — XX54), 

XXig = 2 (XX50 ~ XX51 + XX54). 


(D.X) 


Because of tXie symmetries under the exchange of the momenta qi O —( 72 , only a subset of X2 of these functions 
appears in the master formula (4.28) for {g — 2)^: 


XXi — XXi + (71 • 92XX47, 

1X2 = XX2 — 2 (91 + 9 i ■ 92) (2XX47 — XX50 — XX51 — XX54), 

1X3 = XX4 + {^ q \ + qi ■ 52) XXig + ((7i • (72 + 92) 1X20 

+ (9? + 9 i • 92) (91 • 92 + 92) XX31 — 2 (2XX47 — XX50 — XX51) +2 (91 — 92) XX54, 

1X4 = XX5 — (71 • 92XX21 + 2 (91 ■ 92 + 92) (2XX22 — 2^1 • 92XX33 + XX50 + XX51 — XX54) — 9IXX47, 

1X5 = XX7 — XX19 — ((71 • <72 + 9I) XX31, 
iXe = XXg — XX22 + 9i ■ 92XX33, 

1X7 = XXio — XX21 — (91 • 92 + 92) XX33) 

IXg = XX16 + XX47 + XX54, 
iXg = XX17 + XX47 — XX50 — XX51, 

XilO = 2 (XX39 + XX40 + XX 4 g) , 

fill = XX42 — XX47 + - (XX50 + XX51 + XX54), 

1X12 = 2 (XX50 — XX51 + XX54). 


(D.2) 
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E Integral kernels 

E.l Intermediate kernels 

After calculating the trace and performing the contraction of the Lorentz indices in (4.19), one finds the following 
integral kernels: 


Tiiqi,q2;p) = 


= -3 {{qi ■ q^r - qU) - ^ql^p ■ q,f - • q^^ 

- ^qlp ■ 92 {ql + 9i • 92 ) + P • 9i Q ( 9 ? + 9i • 92 ) 92 + yP • 9291 ' 92^ , 

= ((91 • 92)^ - qlql) ml - ^ql{p-q 2 f +p- 9i (^^P ' 929i ' 92 - ^qjqi ■ 92^ 

+ P- 92 • q2)A , 


^ 2 ( 91 , 92 ;^) = 


^3(91,92;^) = -3 ((91 • 92)^ - 9?9l) ml - - ql { p-qif + -9?4'' 929i ' 92 
+ P • 9i ^“4(7^92 + 3 (91 • 92)^ + 3P • 9291 • 92^ 5 

^4(91,92;^) = ^9i • 92 (9? + 92 + 291 • 92) + ^ql{p ■ qif + ^9?(P • 92)^ 

- • 9iP • 92 (9? + 92 + 491 • 92) , 

'+5(9i,92;p) = ^(P- 9i)^9l - (9? + qi • 52) 9I + ^9?(p ' q^f - ^P • 929i ' 92 (39? + ‘2.qi ■ 92) 

+ P- 9i ^9? (39? + 291 • 92) - ^P- 92 (91 • 92 - 9?)^ , 

76(91,92;p) = ^(p • 92 )^ 9 ? - \ml {ql + gi • ga) 9? - ^P ' 92 (39? + 2gi • ga) 9? 


'75(9i,92;p) = 


r7(9i,92;p) = 


'78(91,92;p) = 


'+ 9 ( 91 ,92;p) = 


T’io(9i,92;p) = 


+ ^92(p • 9i)^ + P • 9i Qp • 92 (9? - 9i • 92) + ^9i • 92 (39? + 2gi • ga)^ , 

= ^9i • 92 ( 29 ? + ( 9 ? + 491 • 92 ) 9? + (91 • 92 )^) - ^9?(P • 9i)^9i • 92 

4 8 

+ 39 ?(p • 92)^ (29? + 9i • 92) - 3P • 9iP • 92 (9? + 391 • 929? + (9i ' 92)^) , 

: ^9i • 92 (29? + (9? + 491 • 92) g? + (gi • q 2 f) ml - ^g?(p • q^fqi ■ q 2 
4 8 

+ 39?(p • 9i)^ (2g? + 9i • 92) - 3P • 9iP • 92 (9? + 39i ' 929? + (gi • 92)^) , 

= - ^ (9? + 9i • 92) (29? + (9? + 491 • 92) 9? + (91 • 92)^) 

+ ^9?(P • 92)^ (9? + 9i • 92) - ^P • 92 (9? + 9i • 92) (9? + 9i • 92) (39? + 2gi • ga) 

+ ^(P • 9i)^ (29? + (g? + 5gi • ga) g? + 2(gi • ga)^) 

+ P • 9i Q (3g? + 2gi • ga) (g? + gi • 92)^ + ^P • 92 (g? (g? + 2gi • ga) - g?9i • 92)^ , 
= - ^ (9? + 9i • 92) (g?9? + (91 • 92)^) ml + ^g?(p • ga)^ (g? - gi • 92) 

- ^9?P • 9291 • 92 (39? + 2gi • ga) + ^(p • 91)^ ((9? + gi • 92) 9? + 2(gi • ga)^) 

+ P • 9i (3g? + 2gi • ga) (gi • ga)^ + ^P • 92 (gi • 92 - 9?) 9i • 92^ , 
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Tii{qi,q2\p) 

Ti2{qi,q2\p) 


Ti3{qi,q2;p) 

Ti4{qi,q2;p) 

Ti5{qi,q2;p) 

TiQ{qi,q2\p) 

Tn{qi,q2\p) 

Ti8{qi,q2\p) 

Ti9{qi,q2\p) 


-- -^mlqi ■ q2 {ql + qi ■ <72) “ 3P ' 92 ((gf + 2gi • 52) ql + 2(gi • 92)^) ql 

+ p-qi Q?’' 92 {ql + qi ■ <72) ql + ^qi ■ q2 {{ql + 2 qi ■ ^2) ql + 2 (gi • 92)^)^ , 


= - 3 (9I + 9i • 92) (2gi + (g| + 4gi • 52) ql + (qi ■ q2f) ml 
+ ■ 9i)^ (9! + 9i • 92) - ^73 • 92 (9? + 9i • 92)^ (39! + 291 • 92) 

+ 3(73 • 92)^ (291 + (92 + 591 • 92) 9i + 2((7 i • 92)^) 

+ 73-91^ (9? + 9i • 92) (9I + 9i • 92) (39! + 291 • 92) 

+ ^73 • 92 (91 + 291 • 929i - 9291 • 92) ^ , 

= - ^ (9I + 9i • 92) (9?9i + (91 • 92)^) ml + ^9|(p • qif {ql - 91 • 92) 

- ^73 • 92(91 • 92)^ (39I + 291 • 92) + ^(73 • 92)^ ((9! + 9i • 92) 9? + 2(91 • 92)^) 

+ ?3 • 9i (^^91 • 92 (39I + 291 • 92) 92 + ^73 • 9291 • 92 (91 -92-92)^ ) 

= ~^mlqi ■ 92 (92 + 9i ■ 92) 92 +73 • 92 ^—2(91 ' 92)^ — ^92 (9i + 29i • 92) 9i ' 92^ 

+ 73-91 Q 73 • 92 ( 9 I + 9i • 92 ) 9 I + ^ (( 9 ? + 291 • 92 ) 9 I + 2(91 • 92 )^) 9i^ , 

= ^ml {ql + 9i • 92) (92 + 9i ' 92) (9? + 92 + 29i ' 92) - ^73'9i 73'92 (9? + 92 + 29i ' 92)" 
= ^ (9? + 92 + 9i • 92) (9?92 - (91 • 92)^) 333.^ - ^92(73 • 9i)^ {ql + 92 + 9i • 92) 

~ ^^ l^P ' '^2)^ (91 + 92 + 9i ■ 92) + -^P ■ qiP ■ 9291 • 92 (9? + 92 + 9i • 92) 5 
= -\^ Ip ■ 92 {ql + 9l • 92) 92 + ^373^ (9^9^ - (91 • 92)^) 92 

- ^(73 • 9i)^ {ql + 6 (9? + 9i • 92)) 9i + ^9?(73 • 92)^ (9I + 291 • 92) 

+ ?3 • 9i ^391 (92 + 9i • 92) 92 + 373 • 9291 • 92 (39 i + 291 -92)^ ) 

= “^9273 • 92 (9? + 9i • 92) 9? + ^373^ {qlql - (91 • 92)^) 9? 

- ^(73 • 92)^ {ql + 6 (9^ + 9i • 92)) 9? + ^92(73 • 9i)^ (9? + 291 • 92) 

+ 73-91 ^91 {ql + 9i - 92) 92 + ^73 - 9291 - 92 (392 + 29i - 92)^ , 

= (9? - 92) (9? (491 - 92 + ql ) + 91-92 (791 - 92 + 492)) 

- ^73 - 92 {ql - ql ) 9 i - 92 (9? + 9 i - 92) + ^92(73 - 9 i)^ {ql + 291 - 92 + ql ) 

- ^ 9 ? (73 - 92)^ (9? + 291 - 92 + 9 i) 

+ ?3 - 9 i (^y 73 - 92 (9? - ql ) {ql + 291 - 92 + 92) + ^ (9? - ql ) 9 i - 92 (91 - 92 + 92)^ ■ 


(E.l) 
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E.2 Kernels for the master formula 


In the master formula (4.28), the following integral kernels appear: 

QIt (af - 1) (erf + 5) + QIt (erf - 1) (af + 5) + AQ 1 Q 2 (erf + af - 2) - 


2QiQ2Qlml 


8 (r2 - 1 ) 4 


Qi (gf - 1 ) (Qit (gf + l) + 4(^2 (r^ - l)) - 4Tm^ 8 (r^ - l) (2w^ - Q|) 

QiQ2Qlml Qlnil 

If 2 (af + af - 2 ) Qit (erf - l) (erf + 7) 8 t 


lis = ^ - 


2Q2ml 


Q 2 T (af - 1) (erf + 7) Qi (1 - erf) Qi (l - erf) 2 2 


2Qito2 
4 8r \ 
m2 QiQ2y ’ 




Qfml 


Q? ^Qi 


^ J_ / 4 (r2 (gf - 1) + o-f - 1) ^ Qi-r (erf - 5) (erf - l) 4r 

Qi y Q2ml QiQ 

_ Q 2 T (erf - 3 ) (erf - 1 ) 2Qi (erf - l) ^ ^ 

Qim2 Qlml Ql 

f_8Q^ 

O O + 16^ Y 


T^ = 7^\QI 


t 2 (erf - 1) (erf + 3) + 4 (erf + erf - 2) 4 \ Qir^ (erf - 5) (erf - l) 


Qf T (af - 1 ) (erf + 5 ) 


Q2to2 


— 4t^ H- X I Ql 


Q 2 T (erf + 5erf — 6) 12t\ 2Qf (erf — l) 


Qim2 


.3 I ^Q^t 


T + T - 


+ Ql 32r2 - 


4Qi (r2 + 1) 


4.n3 _ ^Q± 

^'IQ 2 ml J ml 


Tfi = ^ 


1 f Ql (t 2 ((erf — 22 ) af — 8erf + 29) + 2 (—5erf + erf + 4)) 


Q 2 T (2t'^ (^(erf - 3 )^ - 4erf ^ - 26erf + erf (erf - 12) + 37^ 4 .^ 


Qi (t 2 (—Serf + erf (5erf — 26) + 29) — 4 (erf + 2erf — 3)) Qfr (erf — 9) (af — l) 

O O 


2Q2ml 


QIt {^2 - 9) (^f - 1 ) 8 Q 2 T 2 Qi (1 - erf) ^ 

Qim2 Ql QimS Ql 


Qim2 Ql ( 

Q 2 Q? (8r3 + 22r) Qj {8 t^ - 2) 


Qi ( 8 t 2 + 4 ) 


Qi (8 t3 + 34t) 


Qi (36 t2 + 18) 


- 8 Q 2 T (r^ + 5 ) 


- 8 (r2 + 1) 


- 16<?“ (2r“ + 1) - 
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T. = 


1 f Qi {2 (af + ct|^ — 2) — ((erf + 10) erf + Saf — 19)) 


2ml 


+ Qi 


Q 2 T (2t^ (erf — 5) (af — 1) — 2crf + erf (erf + 4) — 3) 4 t 


2m2 


(o'® - 5) {(^2 - 1) , Qi't (c^f - 9) (c^f - 1) , , 2 


2m2 


2Q2ml 


+ 4t2 




Q 2 Q 1 (8r^ + 6r) 


mi- 


+ Qi 


f 2QlT 

V 


- 8 Q 2 +'r)^ 


+ 


Qi (8t2 - 2) 


+ Ql 


2 Ql (6t2 - 1) 


- 8 (r^ + 1) 


<32 


^8 “ ^2 I Ql I /-)2 


Q 


2, 4 2 (2r^ + 1) (erf + erf - 2) \ , g / 4r 4g2T- (r^ + l) (erf - l) 

\Q2 ml 


^Q2 "V 

6QiT (erf - 1) Qi (2 - 2erf) 
02^2 Qim2 




'Qi (8t2+4) 




g2 (8r2 + 4) 


— 16t' 


Tg = g2 / '"l 1 , ^2 1 


girn2 


g?m2 Q2g2 


+x - 


2Q2 8g2T SQiT 


m,, 


+ 


8(r^ + l)), 


Tio = 


Qi Q 2 

1 ( Qi ('^^ (^■f “ 1) ('^f + 3) + 2 (erf + erf - 2)) g^T (af - 1) (erf + 3) 


g 




gim2 


gi ('e^ (erf - 1) (erf + 3) + 2 (erf + erf - 2)) Q\t (erf - l) (erf + 3) 


mf, 


Q2ml 


+ Qi 


St g2r ((erf + 4) af + af (erf + 4) — 10) \ SQ 2 T 


^Q2 ml j Qi 

+ st ^ + x (-leg? (r^ -1) - i6g2gir (t^ -1) - leg? (r^ -1)) 


+ x 


/4g2giT ^ 4g? ^ igi 

V rn2 


2 

_ _ 2 

ml ml 




1 f g? ("T^ (erf - 5) (erf - 1) - 2 (erf + Serf - 4)) 


Ql ('e^ ((2 - Serf) erf + 1) + erf - Serf + 2) egfr (erf - l) 


mt 


Q2ml 


+ Ql 


+ x 


g2T ((erf — 2) af — 2erf (Serf + 8) + 23) 12r\ Q^t (af — l)" 


2ml 


Q 2 


2Qiml 


' QlQl (St^ +10) 28 Q 2 QIT 12Qi 2Ql 


rat 


rut 


+ -^--^+Ql (S-St^) +8t2 


m 


m 
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Ti2 =- 


QlT ^ - 7) Q2T (erf ^ - 7) 2 QI 


2Q 


where 


4Q2to 2 4Qim2 Qlml Qlnfi^ Q\ Q\ 

\ (Q\ (dr^erf + 2T^crf + Serf + erf) Qf (2T^crf + 4r^CTf + erf + Serf) 


Q3 


2wf. 


llQirerf 

2Q2to2 


2Qiml 


+ 


2mf, 

H' 


X 


Q] 


\ ■ 

- 4Qf (7t^ + 1)+QI (28t^ + 4) + ^ , 



2 ml 


-^ 2 ^) 

2Qfaf 

2Qiaf 


Qlml 

Qlml 

/ 44 Q 2 T 

16r\ 

14Q\Qit 

V 


ml 


X = 


1 


z = 


Q1Q2X 

Q1Q2 


arctan 


zx 


1 — ZT 


= \/r^ 


4to?, 


(1 - eri )(1 - er2 ), a, = \1 + 


4to?, 


Q: 


2 ’ 


Q\ — Q\ 4 " ‘2Q1Q2T + Qf 


(E.2) 


(E.3) 


F Dispersive representations of scalar functions 

F.l Double-spectral representation of basis functions 

In order to derive the double-spectral representation of the basis functions 11^ that is implied by the Mandelstam 
representation of the BTT functions 11^, we have to apply subtractions to the double-spectral integrals. For the 
subtraction of double-spectral integrals at u = uq) we use the following relations: 


If 11 If 

/ ds'du'— -^- Psu{s',u') = ^ / ds'du' 


s' — s W — u 


-t: ( dt'du' 


1 


1 


t' — tu' — u 


Ptu{t',u') = :^ J dt'du' 


1 Psu ^ ) 

^ J S' — S u' — Uq 

M - UO /■ , , , , 1 1 Psu{s',u') 

TT^ J s' — S u' — U u' — Uq 

1 Ptu (t ^ U ^ 


U — Uq 


t' — t u' — Uq 

1 1 


dt du 


/ ds'dt' 


s' — st' — t 


Pst(s',t') =[ ds'dt' — 

TT^ J S 


Ptu (t ^ U ^ 
t' — tu' — U U' — Uq 


1 Pst{s',t') 


4 / ds'dt'- -^- 

— S u' — Uq TT^ j t' — t u' — Uq 


U — Uq 


ds'dt' 


1 1 Pst{s',t') 


7T^ 


s' — st' — t u' — Uq ' (F-1) 

The first two relations appear like subtractions of a single-dispersion integral, while the (st)-double spectral 
integral involves two subtraction terms. Analogous relations are used for subtractions at s = sq or t = to- 

We insert in (5.23) once- or twice-subtracted double-dispersion relations for the scalar functions 11^ to obtain 
double-spectral representations for the basis functions 11^. The functions requiring one subtraction are given by 

1 ^,,DT{t';u,y 

TT 


-2^3 • 54119 = (s - Sb) - / ds 

JAMl 


1 




1 f° 

+ {u — Ub) { - / dt 

JiMl 


s' — s 

,D'^f{t']Sb) 


lAMl 


t' -t 


-I- (s - Sb){u -Ub)\^ 


TT^ 


t' -t 
ds'dt' 


1 f- ^ D2n<sb) 


u — u 

P9;st{s ,t) 1 f P 'p ' P9-,su{s ,U ) 

is'-s)it'-t) J "" “ {s'-s){u'-u) 


dt'du' 


P9;tu{t' ,u') 

[t' — t){u' — u) 
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2qi ■ 921136 = (s - So) i - 

71" ./4MB 


ds'^SM^ + 1 

s' — s 


/'■ J4M2 t' — t j 


- (u-Ub) <1 - I dt 

./'■ J4M2 


, 1)37 (t'; So) 


i'-t 


1 

Tl" .74MB 


du' 


+ (s - Sa){,U - Ub) 


' ds'dt' 


u' — u 

1 


TT^ 


, / , , PZQ-su{s',u') 
{s' — s){t' — t) ' TT^ y {s' — s){u' — u) 




— t){u' — w) J ’ 


-293 -941139 = (5-/^.) I ^ r 

I/'' JAMl 


-(s-5a){- r 

I /'■ J4MB 


i'-t 




1 /•“ 
Tl" AmB 

00 


du' 


,D2i;{u';s,) 


u' — U 




du‘ 


,D-bnu';sb) 


JiMl 


u' — U 


. { \{ ^ ^ U 'UV P39;st{s ,t) 1 f j ij ! P39;su{s ,U ) 

+ (. - .4(. - . 9 1 y * * i^rr^jijTTTj +;? y * 

+ 47 df'iiu' 


{f — t){u' — u) 


-293 • 94n40 = (s - S 6 ) I - r ds' ^50 1 

I/'' JAMZ 




-{t-tb)[- r dt 

I 7’’ d4MB 


5 — S TT J4^M2 U — U 

,D TA{t'-,Sb) ^ 1 r ^^y D3Ain'-,Sb 
t' -t IT J^J ^2 u' - U 


\ I \( 4 - 'U+' PA 0 -,st{s ,t) 1 fill! PA 0 -,su{s , U ) 


+ 4 f dt'du' 


{t' — t){u' — u) J ’ 


(F.2) 


where the discontinuities are defined as (note that in all second integrals, the variable in the denominator should 
be traded for the integration variable by using s' +1' + u' = E) 


Df’''{s';u^) :=- [ 1 f p^■st{s',t') 

TT j U' — Ur, 


TT J U' - Ur 

^ fu fPi',tu{t,u) 1 C fPi-^st{s ,t) 


Dl’''{t';ux) ■= - du' 


u — Ur 


- - I ds' 
TT . 


U — Ur 


Df{s'\tr) := 


1 f ^^,Pi-,st{s',t') 1 f ^^JP^■su{s'■,u') 


t' - tr 


du 


t' - tr 


D^’\u'-,tr) ■■= ~ J 


1 f Pi-,tu{t jU ) ^ f ds' ^ 


t' — tr 


t' — tr 


D^^{t';sr) :=- f _ 1 / du' 

TT J s' - Sr TT J s' - Sr 

TrU-,S / / \ ^ f u fPi',su{s ,u) 1 f /Pi-tu{t,u) 

D, {u; Sr) := - ds - - - dt — - -, 

TT J s' — Sr TT j s' — Sr 


(F.3) 


and the double-spectral densities are given by 


/ ' +'^ P9;st(® )7 ) P22-,st{s ,t ) 

P9-st{s ,t ) := - - -^- - -, 

u' — Ub s' — Sb 

PA3;st{s ^t ) P37\st{s jt ) 


P36;st (Z’ J 7 ) 1 = 
P39;st (Z’ J 7 ) 1 = 
P40;st (Z’ ? 7 ) 1 = 


W — Ub 


PA9;st{^ jt ) P54;si(s ,t ) 


S' - Sb 


P50;st{s' ,t') P5A;st{s',t') 


t' - tb 


s' - Sb 


(F.4) 
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The densities of the (su)- and (tM)-double spectral regions are defined analogously. 
The more complicated functions requiring two subtractions are 



t' -t 




^ f P7;su{s,u) 


H- - I ds dv! 

'TT^ 


{s' - s){u' - u) TT^ J {t'-t){u'-u) 



'19 


t' -t 


dt‘ 




'AMI 


t' -t 




1 fill/ Pl9-,su{s ,U ) 

H- K ds du 

TT^ J 

-4gi • 9293 • 941121 = (s - Sa)(s - Sb) 


du‘ 


1 —, -rw;-r H- 3 - / dt du -r^--r 

(s' — s)(m' — u) TT^ j {f — t){u' — u) 

r 

'4M2 t' — t j 


ds' 


>D^2i{s';Ub) , 1 


4 M2 


r I f°° 

+ {S-Sa){s-Sb){u-Ub)<- ds 

I J4M2 


s' — S TT 

s;u,u/ /. 






+ - 


dt 




TI" J4M2 


t'-t 



1 /",/,/ P21 ;s«('S J M ) 


H-;r / ds du 


{s' — s){u' — u) TT^ 


H-FT / dt du 


/ , / P21-,tu{t',u') 


{t' — t){u' — u) 


(F.5) 


where 


/ / ,/x_ ‘^P7;st{s', t') ^ P9l-st{s', t') 

^ ’ {f - tb){u' - Ua) s' - Sb 

/ / ,/\ 2pi9;si(s ,t ) P31;st(s ) 

Pl9;st(s ,t ) := -r H- - -, 

{U' - Ua){u' - Ub) s' - Sb 

/ / ,/x_ P21;at(s',^') _ P 22 -,st{s', t') 

^ ’ ■ {u' - Ua){u' - Ub) {s' - Sa){s' - Sb)' 


(F.6) 


Again, the densities of the {su)- and (tM)-double spectral regions are defined analogously. The discontinuities 
of the twice-subtracted integrals are: 


Df’''{s';too,Uy) ■■= - I du'-— 


1 


/ Pi\su ('^ J ^ ) 


- - I dt' 


Pi\st ('^ J ^ ) 


TT , 


{P - tx){u' - Uy) TT J {P -tx){u' -Uy)' 


48 




’“(t': 

) ^X-) 

Uy) 


“(s'; 

’^Xl 

Uy) 

DPu, 

“(f'; 

'^x-i 

Uy) 


’*(t': 

i Sx 7 

Sy) 

Du;s, 

S / / 

{u : 

i Sx 7 

Sy) 


- 1 rfn' - 

TT 7 (f - tx)[u' - Uy) 

1 f Pi;su{^ 5) 

IT J {u' - Ux){u' - Uy) 

1- [ , / Pt-tu{'^',U') 

TT J ^ {u' - Ux){u' - Uy) 

1 f , > 

TT J ^ {s' - S^){s' - Sy) 

1 f Pi;su{s jU ) 

TT J {s' - S^){s' - Sy) 


1 [do' 

TT J ^ {t' - tx){u' - Uy)' 

1 [ Pi-,st{s' 

TT J {u' - Ux){u' - Uy)' 

_1 f , ,_PdMs\f)__ 
J '*(«'- Ux){u' - Uy) ' 

1 /■ . , P^■tu^t',u') 

ttJ {s'- s,){s'- Sy)' 

^ f Pi-,tu{i'^u') 

TT y (s' - Sa;)(s' - Sy)' 


The signs are determined by the second subtraction. 


(F.7) 


F.2 Dispersion relation for the pion-box input for (g — 2)^ 

In the limit k —>■ 0, the scalar functions I{i{Qi, Q 2 , r) that are required as an input in the master formula (4.28) 
can be obtained from the basis functions 11^. Due to the presence of kinematic singularities in the basis functions, 
it is important to reach the limit A: —)■ 0 in a very specific way in order to correctly identify the BTT functions 
Ili. One possibility is the following: 


ni((5i, <52, t) 
Tl2{Ql, Q2 ,t) 
^4{Qi, Q 2 , t) 
n5((5i, <52, r) 

^7{QiiQ2,t) 

n9((5i, <52, t) 

nio(( 5 i 7 Q 2 , t) 
ni6((5i7 Q2,t) 
ni7((5i7 <527 t) 
ni9((5i7 Q2,t) 
n2o((5i7 Q 27 t) 
^ 2 i{Qi, Q 2 , t) 
^22{Qi, Q2,t) 

^3i{Qi,Q2,t) 
^33{QiiQ2,t) 
n39((5l7 <52, t) 
n4o(Qi7 Q2tt) 
n42((5i7 <527 t) 
n46((5i7 Q 27 t) 


lim III, 
fc->0 

lim Ci 4 [llil 
fc-^0 '■ ■' 

lim 114 , 
fc-^0 


limCi4 [ 114 ] 7 

lim ( lim IIj j, 
k^0\t^tb J 

lim (Ci3[C23[ lim nr]]) 
fc->o V O-j-tb -'v 

lim fc23[ lim IIt]), 
lini ^C 24 [Ci 4 [ ^lim Er] ] 
lim(c24[Ci3[^limn7]] 


fc-^0 

lim ( lim Ilig ), 

fc—>-0 \ u^Ub 


lim (C 23 
fc->o V 


7 . I 9l • <7293 • 

lim I-1121 


( 


L s^sb V qi ■ 94^2 • 93 

im (C23[ lim ni 9 l ), 

. / / 91-9293 - 9474 

im lim I-1121 

^0 V s^sb \ qi ■ ( 74(72 • 93 

lim f lim ( hJ) 

fc-i-o V s^sb \qi ■ q 4 q 2 -93 / / 

93 • 94 


lim 
k 

lim 
k 


lim (C 23 lim ( 
fc—^0 \ L s^Sb \ 


-n 


Ls^Sb \qi ■ 9492 • 93 
92 


19 


lim fCi 4 lim 
k—bO V L s-ls^ V ( 7 i . 


fc->0 

lim (Ci 
fc-i-o V 


-E; 


9i • 94 


36 


--14 


--24 


lim (Ci 2 
fc-i-o V 

lim (Ci 4 C 23 lim ^^ Eoe) ) 

fc-i-o \ L L s-i-Sa V Oi • (74 / J J / 


9i 

• 92 

9i 

• 94 

9i 

• 92 

9i 

• 94 

9i 

• 92 

9i 

• 94 
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^i7{QliQ2,T) 

n5o((5i7 Q 2 , t) 
n5i((5i, <52, t) 
n54((5i, <52, t) 


lim 1 

{C2 a\ 

<^13 

lim ( 

+1 ■ 92 

fc-^O 

\ 1 


. S-i-Sa \ 

■9i • 94 

lim 1 

{ lim n4n 1 , 


fc-^O 

\ t — 

b 

J 


lim 1 
fc-^O 

(^12 

lim 

t — yt) 

^n4o]) 

5 

lim 1 
fc-^O 

(^23 

C 24 \ 

lim ft 

■ t—Vtb 

»11)- 


(F.8) 


The scalar products in these equations follow from the general relations between the BTT functions and the 
basis functions (B.2). 

Due to the crossing relations, many other equivalent prescriptions exist. The crossing operator Cij is under¬ 
stood to act on the Mandelstam variables s, t, and u as well as on the qf. In particular, the limits Sa, Sb, ■ • ■ 
are affected as well, which ensures that the limits taken in the argument of the crossing operators lead indeed 
to a limit compatible with A: —)■ 0. 

We insert the double-spectral representations for the basis functions (F.2) and (F.5) into these relations. The 
crossing operator acts on the external Mandelstam variable in the Cauchy kernels, but the primed integration 
variables are understood to be unaffected. We perform the action of the crossing operators on the Cauchy 
kernels and rename the integration variables where it is convenient. Finally, we obtain the following dispersion 
relations: 


Q 2 , t) 
n2(Qi, Q 2 , t) 
n4(Qi,Q2,r) 
IIslQi; Q 2 , t) 
n7(Qi, Q 2 , t) 
n9(Qi, Q 2 , t) 
nio(Qi; Q2,t) 
ni6(Qi; Q 2 , t) 

ni7(Qi, Q 2 , t) 
ni9(Qi, Q 2 , t) 
n2o(Qi7 Q 2 , t) 
n2i(Qi, Q 2 , t) 
n22(Qi, Q 2 , t) 


lim 

lim 

lim 

lim 

<?|^o 

lim 

lim 

<?|^o 

lim 

lim 

lim 

lim 

lim 

lim 

lim 


JiMl 

'1 r 
'1 r 
'1 r 

'1 r°° 

/*oo 


wS . 2 'I 

^ J4M^ I' ”1" W 2 


ds 


ds 


ds 


s' + Qg VI J 4 M 2 ^ w ^;2 y 

,Ci 4[D*h"(5>b)] , 1 r , Cu[Dr{t';ub)] \ 
s' + Ql t' + Ql )' 

, Dris'lu,) 1 /■“ Dr{t';u,) \ 

' s' + Ql t' + Ql )' 

, c^4DTis';u,)] 1 Cu[Drit';u4] \ 

s' + Ql ^Jami t' + Ql 


1 f- Driu';tQ \ 
s' + Ql u' + Ql )’ 

j Ci3[C23[Df{s';h)]] ^^, C33[C23[Df{t';h)]] 

s' + Ql TT J 4 M 2 t' + Ql 


1 

Jam^ 

'1 r 

Jam^ 

'1 r 

^ Jam^ 

'1 r 

^ Jam^ 
Jam^ 

, ^ Jami 


s' TT 74M2 i' + Q 2 

., C23[Df{t';h)] 1 , C23[Dfiu';h)] \ 

t' + Ql TT i4M2 "" u' + Ql )' 

r . , C2A\CAA[DT{s'-h)]] , 1 r ,., C2A\ClA[Df{t';tQ]] 
JAMl s' + Q3 ^ JaMZ t' + Ql 


Jam^ 

vC24[Cl3[Df (s';4)]] ,1 r . ,C2A[Cl3[Df{u'-h)]] 

s' + Ql ^Jami "" u' + Q\ 

, D^4{s'-,uQ ^ 1 /■“ 

s' + Ql 


ds' 

ds' 

dt 


' + QI 


-Ql 


Di 


t' + Ql 


'■ JAMl 

Q 2 y 

1 

-k - 

f>00 

C23[D^22\u'-,Sb)] 

TT , 

Iam^ 

u' + Ql 

1 

■ -k - 

nOD 

C23[D^4{t';uQ] 

TT 

Jam^ 

t' + Ql 

nOO 

. v^22'K;Sb)\ 

Jami 


u' + Ql J 
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n3i(Qi,Q2,T) = lim 


1 


dt 


\ TT J4M2 

Y roo 


^ 1 /■“ 

t' + Qi Am 2 


du‘ 


' + Q? 


n33(Qi,Q2,r) = lim - 


fis' 


q|->0 y TT J4M2 
2 /•CO 


n39(Qi, Q 2 , t) = lim - 


ds' 


(Jl-s-O \ TT J 4 M^ 
1 


C23[D^3'lis';Sb)] 

+ - 1 

^00 

C23mnn';sb)] 

s' + 

77 J 

4M2 

u' + Ql 

Ci4D*^^{s';Sa)] 

1 

poo 

^ ,Cu[D-i^{n'-,Sa)] 

s' + (3§ 

! 

4 M2 

u' + Qi 


n 4 o(Qi,Q 2 ,r) = lim - 


ds 


9|->0 y TT J4JV42 


,Cl2 [Cl4[J37(g';ga)]] 
s' + Ql 


/ 1 

n 42 (Ql,Q 2 ,T) = lim - 


^^/^2[C24[i^3T(5';Sa) 


1 

TT , 

1 


dt‘ 


,Ci2[Ci4[i^3T(i';Sa)]] 


'4M2 


t' + (32 


,Ci 2 [C 2 a[Dz 7 ®o)]] 


<?4-s-0 y TT J4M2 s' + <33 

,i;s / 


du‘ 




u' + Q 


1 


n46(<3i,<32,r) = lim - 

9|^r - 

n47(<3l,<32,T) = lim 


dt 


j Ci4 [C23 [-D37 (t'; Sq 


(jI-S-O y TT J 4 JV 42 
1 

\^7r ^4^2 

n5o(<3i,<32,'r) = lim ( - [ 

ql-^0 y TT d4JV42 


t' + <3 


“ C24pl3[i?|^(i^ 

t' + Ql 


1 

d4MJ 

1 


du‘ 


,Cu[C 23 m\u';Sa)]] 


du‘ 


77 


/1 

n5i(<3i,<32,r) = hm - 


j / ^50 ('^'5 ^fe) 

s' + Qi 


1 r°° 

^ Am? 


'4M? 


w' + (3i 
w' + (3i 


,C24[Ci3[iAfVlM] 

/ _L /02 




ti' 




(jI-s-O y 77 J4M? 
1 


Qi 


1 

77 


n54(<3i,<32,r) = lim - 


dt 


,C 23 [C 24 [D'^o\t';tb)]] 


(Jl-S-O y 77 J4JV42 


t' + <3 


/ dt 

'4M? 

- / 

77 J4M? 


Q? y’ 

,Ci2[i?^6‘A;4)]' 


du‘ 


t' Q2 J 

,C23p24[£’50 K;^b)]] 


Q? 
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The limit fc —)■ 0 is already taken in the Cauchy kernels. As the numerators only depend on the integration 
variables and the qf, the remaining limit to be taken is <74 —)■ 0. 

For Hi, 112, 114, and IIs we have randomly chosen a fixed-u dispersion relation. Because the first six basis 
functions are not affected by kinematic singularities, we could equally well choose a fixed-s or fixed-t dispersion 
relation or use directly the Mandelstam representation. 


F.3 Kinematic zeros due to crossing antisymmetry in HLbL 

Here, we discuss in more detail the crossing property of the BTT function 049. According to (3.15), n49 and 
its crossed versions are the only scalar functions with odd intrinsic crossing properties. Crossing antisymmetry 
implies an additional kinematic zero, but, so far, we have ignored this and written down dispersion relations 
directly for n49. As neither the pion-pole nor the pion-box topologies contribute to n49, the treatment of the 
kinematic zero in n49 can be of relevance only for contributions of higher intermediate states than discussed 
here. We show here in general that the dispersive representation of the basis functions is not affected by the 
presence of the kinematic zero. In particular, the kinematic zero has no influence on the contribution to {g — 2 )^. 

Let us denote 


n49^^(s, t, u) ■— n49(s, t, u] q\, (72, (73, q\)- (F-10) 

The crossing properties (3.15) of 049 are then 

nif'‘(s, t, u) = -n^^3'‘(s, u, t) = -n^^'^3(s, u, t) = n^9'‘^(s, t, u). (f.ii) 

We build symmetric and antisymmetric combinations under exchange of the Mandelstam variables t and u, but 
fixed virtualities: 

5']^^‘‘(s, t, u) := i ^(s, t, u) + '‘(s, u, t)) , 

Alf Aa t, u) := i {Ulf\s, t, u) - Ull^^s, u, t)) . (F.12) 
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They fulfill the crossing relations 


t, u) = u, t) = t, u) = Alf^^{s, t, u) 

and hence exhibit the following kinematic zeros: 

Sll^‘^{s,t,u) = {ql - ql){ql - ql)Slf‘^{sA,u), 

t, u) = {t- ■u)A49(s, t, u). 


(F. 13 ) 


(F. 14 ) 


The functions 5*49 and A49 are symmetric under the individual exchanges of t o u, and g| g|. The 

question arises now if it makes a difference if we write down dispersion relations for 1149 or rather for S'49 and 
A49. We show in the following that the Mandelstam representations are not equal but that the basis functions 
Hi remain unaffected. Therefore, the kinematic zero has no influence on any physical quantity. 

Let us start with the Mandelstam representation of 1149: 




ds' / dt' 


1 


H —K dt' dvL 


TT^ 


(s' — s){t' — t) if- 
{f — t){u' — u)' 


J ^ J ^ (s' - s)('u' - u) 


(F. 15 ) 


The implicit integration limits are the borders of the double-spectral regions. The crossing properties of 1149 
imply the following relations for the double-spectral densities: 


= -plltLis,t) = Pllt!tis,t), 
plZiit,u) = -pl^l{u,t) = -pZH^a) = pZiit,u). 

With the Mandelstam representation for 1149, we obtain for the functions defined in (F. 12 ): 


(F. 16 ) 


rfl234 

^49 


is,t,u) := - 


1 / 1 


-2 

TT^ 


ds' / dt 


dt / du 


:= ^ J ds' J dt 


(s' — s)(t' — t) TT^ 

{t' — t){u' — u) J’ 
,plif4(s',t')-pii?i(s',F) , 1 


-2 

TT^ 


dt' / du' 


{s' — s){t' — t) 7 T- 

{t' — t){u' — u) j 


ds' / du‘ 


,pII %Us',u') + pII%Us',u') 

{s' — s){u' — u) 


H —K ds' / du‘ 


,pll%is',u')-pll%{s',u') 


(s' — s){u' — u) 


(F. 17 ) 


The crossing properties of the double-spectral densities allow us to factor out kinematic zeros. In the case of 
S'49, the combinations of double-spectral densities have kinematic zeros of the form {qf — g|)(g| — gl), which 
can be taken immediately out of the dispersive integrals, hence we find indeed 




TT^ J J (s'-s)(<'-t) 7r2 

7 r 2 y J ^ {f - t){u' -u)]' 


ds' / du‘ 


(s' — s){u' — u) 


(F. 18 ) 


where the (J49 denote the double-spectral densities of S49. 

In the case of 4I49, some algebra is necessary. We subtract in each term the dispersion integrals so that 
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factors of {t — u) multiply the double-spectral contribution: 


TT^ 7 y [s' - s){t' - t) 


= 1 fds'^—- - P 49 %is',t') 

TTj^ s' - STtJ f - 

1 f 1 1 f ^j P49%is'^ t') - piifKg'. t') 

TT J t' -tTT J ^ 2{t' - 

1 Pll%{s',t')-pll%i{s',t') 


- s){t' -t) 2 {P - 


f ds' I du' 


, pII% 7 ',u')-pII%{s',u') 

{s' — s)(u' — u) 

= 1 fds'^—- [ du ' ~ 

ttJ^ s'- sttJ^ u' - 

H— du — - ds -p- 

TT J u'-uttJ 2 (u'- 7 ^) 


\^ [ U ' [ U ' 1 P49?i(s',w') -P49?s«(s'>w') 

TT^ J J {s' - s){u' - u) 2 (u'-77-) 


— f dt'- f d ' '^ 49 ftt(^'. ^0 - pjjfttK: t') 
J J {t' — t){u' — u) 


1 fdt'^^ f^^> Pl^tiit',u')-pl^ 7 {u',t') 
TT J t' — tTT J U' — t' 

1 1 1 f ,uP^tit\u')-p]^l{u',t') 


H— / du 


TT . 


/ -— dt‘ 

W — UTT J t' — U' 

,1234 C+/ ^1234 j./\ 


,11 1 /■, ,riSi( 4 ',«')-(>Si(»'.(') 


(F. 19 ) 


We see that in the sum of these three terms, the one-dimensional s-channel integrals cancel but that the other 
subtraction terms do not. Hence we find: 


4 1234 
"^49 


(»,<,») = (t-«)(l i /‘‘"'/‘‘“F 


, algfKs'.w') 


— s){u' — u) 


+ F fdi' l,n' ««'.»') 


{f — t){u' — u) 


AgW Al-M 


f — t TT , 


W — U 


(F. 20 ) 


where the 0:49 denote the double-spectral densities of 4I49 and where 


4 12347.11 1 [ , 1 [ , ,P^^L{t',u')-p^l{u',t') 

A49 (t):--jds 2 t' + s'-i: + ^■ 


(F. 21 ) 


We learn that writing a Mandelstam representation for n49 is not equivalent to taking explicitly into account 
the kinematic zeros and writing Mandelstam representations for 549 and H49. The difference consists of a 
contribution with t- and u-channel cuts but no double-spectral regions. 

However, due to the redundancy in the set of BTT functions, n49 is not a physical quantity: the HLbL 
tensor is defined by the basis functions Hi, and the dispersive representation for these does not change whether 
one starts from an unsubtracted Mandelstam representation of n49 or £'49/449. To prove this, we need to study 
the impact of the difference terms in (F. 20 ) on the basis functions H^. The double-spectral representations for 
the basis functions involving n4g are given in (F. 2 ). The difference terms in (F. 20 ) obviously result in a shift 
of the single-variable dispersion integrals, while the double-spectral contributions are unchanged. Note that 
the single-variable dispersion integrals involve only fixed-s representations of n49 and crossed versions thereof. 
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They have exactly the form of the difference terms in (F. 20 ), hence the difference between the choices of 1149 or 
S'4g/^49 as the functions to fulfill an unsubtracted Mandelstam representation can be reabsorbed in a shift of 
£>49 and -D49®. This means that the basis functions Hi and the HLbL tensor remain unaffected by the kinematic 
zero in 1149. In particular, the form of the dispersion relation for {g — 2 )^ in (F. 9 ) is unchanged: the functions 
IIso, Hsi, and 1154 are defined by crossed versions of a fixed-s dispersion relation for 1149, where the effect of 
the kinematic zero is reabsorbed in the discontinuities of the t- and u-channel integrals. Since in the end the 
discontinuities will be calculated directly from the unitarity relation, the observation that the kinematic zero in 
1149 does not alter the functional form of the dispersion relation completes the proof that there are no practical 
consequences, neither in the construction of the dispersion relation nor in [g — 2)^. 


G Dispersive representation of loop functions 

G.l Scalar two-point function 



m2 


Figure 9: Bulb diagram 

The scalar two-point function [ 64 ] corresponding to the loop diagram in Fig. 9 is given by 

d'^q 1 




( 27 r )4 [g2 _ ml] [{q + p)2 - ml] 


(G.l) 


We define s := p^. The two-point function has a normal threshold at s = (mi -l-mg)^. According to Cutkosky’s 
rules [ 63 ], the discontinuity is given by: 


2 iAsBo{s) := Bo{s + it) — Bq{s — it) = 


i A^/^(s,mi,m2) 


Stt 


Bq is divergent and satisfies a once-subtracted dispersion relation: 

, AsBo{s') 


Bois) - Bo{ 0 ) = - f 

J (1 

G.2 Scalar three-point function 


ds' 


{mi+m 2 )^ 


{s' — s — it)s' 


(G.2) 


(G.3) 



Gonsider the scalar three-point function [ 64 ] corresponding to the loop diagram in Fig. 10 : 


C'o(£i I £21 £31 



d^q _1_ 

(27r)4 [g2 _ ml][{q + pi)"^ - ml][{q + pi + P2)‘^ - ”^- 3 ]' 


(G.4) 
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We define s := (pi +^2)^ = pi- The scalar three-point function has a normal threshold at s = (mi +7713)^. The 
discontinuity along the corresponding branch cut can again be found with Cutkosky’s rules: 


2*A,Co(s) 


(70(5 -I- ie) — Co(s — ie) 


i mf, m3) 

327 r 2 s 




(q + Pi)'^ - ml' 


Performing the phase space integral leads to: 

. „ , ^ ^ 1 1 1 / 2s(Ai2 +Pi) - (s + Ai3)(s -I- A12) -I- Ai3^(s)Ai2^(s) \ 

W-Ww A2»(A.2+p;)-(.+A,3)(,+A„)-A;f(.)A;f(.))’ 


where A^fc := pf - pi, Aik ■= ml - ml, Aifc(s) := X{s,pl,pl), and Xik{s) := X{s,ml,ml). 
Co satisfies an unsubtracted dispersion relation: 


Co(s) 


1 r 

^ J (mi-t-ms)^ S S ic 


(G. 5 ) 


(G. 6 ) 


(G. 7 ) 


As the anomalous threshold appears this dispersive representation needs to be amended by adding a dispersive 
integral over the anomalous cut [ 55 - 57 ]. 


G.3 Scalar four-point function 


Pi Pi 



Figure 11: Box diagram 

Gonsider the scalar four-point function [ 64 ] corresponding to the diagram in Fig. 11 : 

Do{pj,pl,pl,pl,{Pi + P2f ,{P2 + p-if ,ml,ml,ml,ml) = 

1 f 1 

i J ( 27 r )4 


(G. 8 ) 


We define 


s := (pi -kp2)^ 
t := (p2 +P3)^ 
u := {pi +P3f, 

P :=pi+P2 = -p 3 -p 4 . (G. 9 ) 

In the s-channel, the diagram has a normal threshold at s = (mi -I- m3)^, in the f-channel at t = (m2 -I- m4)^. 
For a fixed value of t, we find the discontinuity along the s-channel cut: 


2zA*I?o(s, t) := Do{s + ie, t) — Do{s — ie, t) 

i X^/^(s,ml,ml) f „ 1 _ 1 

327r2 s J {q + piY - ml {q - Pa)"^ - m\' 


(G.IO) 
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The phase-space integral can be transformed into a dispersive integral, see [ 33 , 76 ]: 


where 


and 


1 1 , 1 

= — / dt - , 

:= aia 2 ± ]/(af - 1) (a^ - 1), 

5(2?7l2 — + Ai2^13 

_ s(2m| - S") -I- - A 34 A 13 

E :=pf +P2 +pI +pI 
E' :=p2 +p2 

Y^// 2 I 2 I 2 I 2 

E :=Ps+P4 + 771^+7713, 

A*fc :=pI-pI, 

A,fc := TO- - TOfc, 

A*fe(s) := \{s,pI,pI), 

\k{s) := X{s,m1,ml). 


(G.ll) 


(G. 12 ) 


(G. 13 ) 


Since Dq satisfies an unsubtracted dispersion relation, it can be written as 

Po{s',t') 


Dois,t) = \r ds' f 


dt' 


(mi+ms)^ Jt+{s') ^ ic){t t ic) 


(G. 14 ) 


where 


Pois,t) = 


8Aif (s) ^yit-t+(s)){t-t-{s)) 


(G. 15 ) 


The validity of the double-spectral representation of the box diagram breaks down with the appearance of 
anomalous thresholds. The possibility to amend it by adding the contribution of a dispersive integral over 
anomalous cuts has not been investigated here, as it is irrelevant for the {g — 2)^ calculation. 

If we do not take a fixed-t dispersion relation as the starting point but a fixed-u dispersion relation, we have 
to take into account that a line of fixed u in the Mandelstam plane will encounter a right- as well as a left-hand 
cut: 


1 /■“ 

Dois,t) = - ds‘ 

J-00 


,a^Do{s',j:-u-s') 


s' — s — le 


where 


A^Do{s',T,-u- s') 

, 2sDo{s'+ ie,J: - u-s') - Do{s'- ie,T, - u-s') 

= 0 {s - (toi -k TO3) )-—- 

2 i 

+ (!(E(m, + S-I-«) 


2 i 


hence 


1 

Dois,t) = - ds‘ 


,A^Do{s',X:-u - s') 1 


s' — s — le 


^ J (m 2 +m 4)2 


dt 


,A^Do{Y.-u-t',t') 
t' — t — ie 


(G. 16 ) 


(G. 17 ) 


(G. 18 ) 
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Cutkosky’s rules lead to the following discontinuities: 


A'^Dois,J:-u - s) = 




1 1 


dt'- 


1 


Stt 7i+(s) {f - Ys + u + s - ie)y^{t' - t+){t' - t“) ’ 

1 


'13 

1 1 


ds' 


Js+it) {s' -Y + u + t- ie)y^{s' - s+)(s' - s”)’ 


(G.19) 


where s^{t) are defined in analogy to t^{s) with the proper permutations of the momenta and masses. 
Now, we note that 


/ dt' / ds' = 


/ ds' / dt' 


and find 


where 


2 pOO pOO 

Do{s,t) = ^ / ds' / dt'X{s,t, s',t'), 


X{s, t, s', t') = 


8AJ^^(s') s ' - s-ie t' -t-ie y^(t' - t+){t' - f”) s' - s + t' -t-ie 


m^{t') V («' - s+{t')){s' - s-{t')) j ■ 


With some (computer) algebra, we can check that 

rl/2 


SO that we find again 


i {t'-t+{s')){t'-t-{s')) 
Xli\t')^ i^'-'^+it')){s'-s-it')) ’ 


Y/ i ' _ _ Pojs 1 O _ 

^ ^ “ (s' - s - ie){t' - t - ie)' 


(G.20) 


(G.21) 


(G.22) 


(G.23) 


(G.24) 
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